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Abstract. This paper presents a geometric description of Lagrangian and Hamiltonian sys- 
tems on Lie affgebroids subject to affine nonholonomic constraints. We define the notion of 
nonholonomically constrained system, and characterize regularity conditions that guarantee 
that the dynamics of the system can be obtained as a suitable projection of the uncon- 
strained dynamics. It is shown that one can define an almost aff-Poisson bracket on the 
constraint AV-bundle, which plays a prominent role in the description of nonholonomic dy- 
namics. Moreover, these developments give a general description of nonholonomic systems 
and the unified treatment permits to study nonholonomic systems after or before reduction 
in the same framework. Also, it is not necessary to distinguish between linear or affine 
constraints and the methods are valid for explicitly time-dependent systems. 
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I. Introduction 

Nonholonomic constraints is one of the more fascinating and applied topics actually. First, 
there are many open questions related with this subject: characterization of the integrability, 
construction of geometric integrators, stabilization, controllability... But morever, there is a 
wide range of applications of this kind of systems in engineering, robotics... (see [3IE1HH] an d 
references therein). 

During the last years, many authors have studied in detail the geometry of nonholonomic 
systems. Some of them developed a geometric formalism for the most typical nonholonomic 
systems, the ones determined by a mechanical Lagrangian, that is, 

L ( V q) = \9{Vq,V q ) - V(q), V q G T q Q , 

where V : Q — > R is the potential energy on the configuration space Q, g is a Riemannian 
metric on Q and, additionally, the system is subjected to linear constraints on the velocities, 
expressed as a nonintegrable distribution on Q [3 [T§] . The usual formalism for these systems 
was the use of adequate projections of the Levi-Civita connection associated to g to obtain the 
equations of motion of the system. Other authors preferred to work on the tangent bundle of the 
configuration space, which permits to introduce more general Lagrangians and different type of 
constraints (linear or nonlinear) . Usually these systems are specified by a Lagrangian function 
L : TQ — > M and a constraint submanifold M of TQ [T3 EH]. Moreover, also a description 
in terms of a nonholonomic bracket was introduced [3 I13[ 136] with considerable applications 
to reduction. In this sense, reduction of nonholonomic systems was intensively studied vi- 
sing different geometric techniques. For instance, introducing modified Ereshmann connections 
(nonholonomic connection) |3j. using a symplectic distribution on the constraint submanifold 
[TJ [35] and distinguishing the different casuistic depending on the 'position' of the Lie group of 
symmetries acting on the system and the nonholonomic distribution or in a Poisson context |22j 
(see [6] and references therein). Moreover, using jet bundle techniques, many authors studied 
the geometry of nonholonomic systems admitting an extension to explicitly time-dependent 
nonholonomic systems [3 U3 [TBI 1301 [33 [31] and ready for extension to nonholonomic field 
theories [57] . 

Today we find an scenario with a very rich theory but with an important lack: a general 
framework unifying the different casuistic (unreduced and reduced equations, systems subjected 
to linear or affine constraints, time-dependent or time-independent systems...). 
The aim of the present paper is cover this lack and present a geometrical framework covering the 
different cases. Obviously, in order to reach this objective it will be necessary to use some new 
and sophisticated techniques, in particular, the concept of a Lie affgebroid (an affine version of 
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a Lie algebroid) and appropriate versions of Lagrangian and nonholonomic mechanics in this 
setting. 

The previous and other motivations coming from other fields (topology, algebraic geometry...) 
have recently caused a lot of interest in the study of Lie algebroids [3T] , which in our setting can 
be thought of as 'generalized tangent bundles', since it allows to consider, in a unified formalism, 
mechanical systems on Lie algebras, integrable distributions, tangent bundles, quotients of 
tangent bundles by Lie groups when an action is considered... In this sense, Lie algebroids can be 
used to give unified geometric descriptions of Hamiltonian and Lagrangian Mechanics. In [38 , 
A. Weinstein introduced Lagrangian systems on a Lie algebroid E — > M by means of the linear 
Poisson structure on the dual bundle E* and a Legendre-type map from E to E* . Motivated by 
this description, E. Martinez [23) developed a geometric formalism on Lie algebroids extending 
the Klein's formalism in ordinary Lagrangian Mechanics on tangent bundles. This line of 
research has been followed in [I7J [24]. In fact, a geometric description of Lagrangian and 
Hamiltonian dynamics on Lie algebroids, in terms of Lagrangian submanifolds of symplectic 
Lie algebroids, was given in [17]. More recently, in [7J (see also [H [57J 28 ) a comprehensive 
treatment of Lagrangian systems on Lie algebroids subject to linear nonholonomic constraints 
was developed. The proposed formalism allows us to treat in a unified way a variety of situations 
for time-independent Lagrangian systems subject to linear nonholonomic constraints (systems 
with symmetry, nonholonomic LL systems, nonholonomic LR systems,...). 
On the other hand, in [5S] an affine version of the notion of a Lie algebroid structure was 
introduced. The resultant geometric object is called a Lie affgebroid structure. Lie affgebroid 
structures may be used to develop a time-dependent version of Lagrange and Hamilton equations 
on Lie algebroids (see [lOj [12l [25] [26] [31]). In fact, one may obtain Lagrange and Hamilton 
equations on Lie affgebroids using a cosymplectic formalism [T2l |2"5] . In the same setting of 
Lie affgebroids, one also may obtain the Hamilton equations using the notion of an aff-Poisson 
structure on an AV-bundle (see [HUH!])- An AV-bundle is an affine bundle of rank 1 modelled 
on the trivial vector bundle and an aff-Poisson bracket on an AV-bundle may be considered as 
an affine version of a Poisson bracket on a manifold (see [9] ) . 

The aim of this paper is develop a geometric description of Lagrangian systems subject to affine 
constraints using the Lie affgebroid theory. The general geometric framework proposed covers 
the most interesting previous methods in the literature. 

The paper is organized as follows. In Section II, we recall several constructions (which will be 
useful in the sequel) about the geometric description of Lagrangian and Hamiltonian Mechanics 
on Lie affgebroids and the equivalence between both formalisms in the hyperregular case. In 
Section III, we introduce the Lagrange-d'Alembert equations for an affine nonholonomic La- 
grangian system on a Lie affgebroid. In Section IV, we discuss the existence and uniqueness 
of solutions for this type of systems. Moreover, we prove that in the regular case the nonholo- 
nomic dynamics can be obtained by different projections from the unconstrained dynamics. In 
Section V, we develop the Hamiltonian description and we discuss the equivalence between the 
Lagrangian and Hamiltonian formalism. We also introduce the nonholonomic bracket which 
gives the evolution of an observable for the nonholonomic dynamics. The nonholonomic bracket 
is an almost aff-Poisson bracket (an aff-Poisson bracket which doesn't satisfy, in general, "the 
Jacobi identity") on the constraint AV-bundle. In Section VI, we apply the results obtained in 
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the paper to two particular cases: a linear nonholonomic Lagrangian system on a Lie algcbroid 
and an standard affine nonholonomic Lagrangian system on the 1-jet bundle of local sections of 
a fibration r : M — ► R. As a consequence, we directly deduce some results obtained in [7l[T4],[T6] 
(see also [1]). In addition, we also analyze (in the Lie affgebroid setting) the reduced equations 
for a well-known example of an affine nonholonomic mechanical system: a homogeneous rolling 
ball without sliding on a rotating table with time-dependent angular velocity. For this purpose, 
we will use a particular class of Lie affgebroids, namely, Atiyah affgebroids. 
The paper ends with our conclusions and a description of future research directions. 

II. LAGRANGIAN AND HAMILTONIAN FORMALISM ON LlE AFFGEBROIDS 

In this section we recall some well-known facts concerning the geometry of Lie algebroids and 
Lie affgebroids. 

II. 1. Lie algebroids. Let E be a vector bundle of rank n over the manifold M of dimension 
m and te '■ E — > M be the vector bundle projection. Denote by T(te) the C°° (M)-module 
of sections of te : E — > M. A Lie algebroid structure ([•, -J E , Pe) on E is a Lie bracket [•, -Je 
on the space T(t~e) and a bundle map pe '■ E — > TM, called the anchor map, such that if we 
also denote by pe ■ T(te) — > X(M) the homomorphism of C°° (M)-modules induced by the 
anchor map then [X, fYj E = f{X, Y] E + p E {X){f)Y, for X, Y G T(t e ) and / G C°°{M). The 
triple (E,\-,-\e,Pe) is called a Lie algebroid over M (see [H]). In such a case, the anchor 
map pe ■ F(t e ) — > X(M) is a homomorphism between the Lie algebras (T(te), [•, -Je) and 
(X(M), [-,.]). 

If (E, \-,-\e,Pe) is a Lie algebroid, one may define a cohomology operator which is called the 
differential of E, d E : T(A k r E ) — > F(A fc+1 r*), as follows 

k 

(d E p)(X Q ,...,X k ) - ^(-ljV^XMAo,...,^,...,^)) 

^ _ _ (1) 
+ X](-l) ,+ V([X i) Xj\ E , x , V V ; X k ), 

i<j 

for fj, 6 F(A fe r*) and X ,...,X k G T(t e ). Moreover, if X G T(t e ), one may introduce, in a 
natural way, the Lie derivative with respect to X, as the operator L\ : T(/\ k T E ) — ► r(A fc T^) 
given by £^ = ix ° d E + d E o i x . 

If E is the standard Lie algebroid TM then the differential d E = d™ is the usual exterior 
differential associated with M. On the other hand, if E is a real Lie algebra g of finite dimension 
then g is a Lie algebroid over a single point and the differential d 3 is the algebraic differential 
of the Lie algebra. 

Now, suppose that (E, [•, -]e, Pe) and (£", [■, -Je', Pe 1 ) are Lie algebroids over M and M', 
respectively, and that F : E — > E' is a vector bundle morphism over the map / : M — > M' . 
Then (F, /) is said to be a Lie algebroid morphism if 

d E {{F,f )*(/)') = (FJ)*(d E ' '</>'), for 0' G r(A fe (r B 0*) and for all fe. 

Note that (F, f )*</>' is the section of the vector bundle A k E* — * M defined by 

(OF, /)*0 , ).(ai, • • ■ , a fc ) = ^ (x) (F( ai ), . . . , F(a fe )), 
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for x G M and ai,...,afc G fr. If (F, /) is a Lie algebroid morphism, / is an injective 
immersion and F\ Ex : E x — > E'^,-. is injective, for all x G Af then (E, [•, -]b,p_e) is said to be 
a Lie subalgebroid of (£", [•, -],e< , p E >). 

If we take local coordinates (x l ) on M and a local basis {e a } of sections of E, then we have the 
corresponding local coordinates (x z ,y a ) on E, where y Q (a) is the a-th coordinate of a G f in 
the given basis. Such coordinates determine local functions p^, C^o on M which contain the 
local information of the Lie algebroid structure, and accordingly they are called the structure 
functions of the Lie algebroid. They are given by 

d 

PE{e a ) = Pa-g-[ and [e a ,ep\ E = C2 /3 e 7 , 
and they satisfy the following relations 



d Pfi J dpi 



^di-^di =p ^ a,,<! 

If / G C°°(M), we have that 



dxi Hf3 dxi ^ af3 * ^ 

cyclic(a,/3,7) 



' dx % /37°"M 



1 dx iPa ' 

where {e a } is the dual basis of {e a }. Moreover, if 9 G r(r^) and — 6* 7 e 7 it follows that 

II. 2. The prolongation of a Lie algebroid over a fibration. In this section, we will recall 
the definition of the Lie algebroid structure on the prolongation of a Lie algebroid over a 
fibration (see [HI HZ]). 

Let (E, [•, -]e,Pe) be a Lie algebroid of rank n over a manifold M of dimension m with vector 
bundle projection t e ■ E — > M and it : Af — » M be a fibration. 

We consider the subset T E M' of E x TAf and the map r| : T E M' -> Af defined by 

T E M' = {(b,t/) G £ x TM'/p E (b) = (Ttt)(v% T E (b,v') = 7r^(«'). 

where T7r : TM' — ► TM is the tangent map to 7r and 7Tm' : TM' — > A/' is the canonical 
projection. Then, t|J : T E M' — > Af is a vector bundle over Af of rank n + dim Af — m which 
admits a Lie algebroid structure ([•, p E ) characterized by 

l(X o 7T, I/'), O 7T, = ° 7T, [C/', V]), pl(X o 7T, C/') = U', 

for all 1,7 6 r(r^) and U', V vector fields which are 7r-projectable to p E (X) and p E (Y), res- 
pectively. (T E M\ [•, ■J E ,p E ) is called i/ie prolongation of the Lie algebroid E over the fibration 
7r or the E-tangent bundle to Af (for more details, see [TTj H"7]). 

An element Z G T E M' is said to be vertical if its projection onto the first factor is zero. 
Therefore, it is of the form (0, v x i), with v x i a tangent vector to Af at x' which is ir- vertical. 
Next, we consider a particular case of the above construction. Let £ be a Lie algebroid over a 
manifold M with vector bundle projection r E : E — > M and T E E* be the prolongation of E over 
the projection r E : E* — > Af T E E* is a Lie algebroid over f * and we can define a canonical 
section X E of the vector bundle (T E E*)* -> f * as follows. If a* G f * and (b, w) G T E E* then 

A B (a*)(M) = a*(b). (2) 
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As is called the Liouville section associated with the Lie algebroid E. 

Now, one may consider the nondegenerate 2-section fig = — d r E Xe of (T E E*)* — ► E*. It is 
clear that d T E He = 0. In other words, Qe is a symplectic section. He is called the canonical 
symplectic section associated with the Lie algebroid E. Using the symplectic section He one 
may introduce a linear Poisson structure II^;. on E*, with linear Poisson bracket {•, - }e* given 
by 

{F, G}e* = tl E (X% E ,x2 E ), for F,G £ C°° (E*), (3) 
where X^ E and X% E are the Hamiltonian sections associated with F and G, that is, i x u E Q,E = 

F 

d rEE 'F and i v n E Vt E = d? E E * G. 

G 

Suppose that (x l ) are local coordinates on an open subset U of M and that {e a } is a local 
basis of sections of the vector bundle t^}{U) — > U as above. Then, {e Q , e a } is a local basis of 
sections of the vector bundle (r^)- 1 ^!,) -1 ^)) -> (r^)- 1 ^), where r^f : T^* -» E* is 
the vector bundle projection and 

?.(a-)=(e (7i(a-)),^ | J, = (o,^ J. 

Here, (x l ,y a ) are the local coordinates on induced by the local coordinates (x 1 ) and the 
dual basis {e a } of {e Q }. Moreover, we have that 

[e Q .,e |3 ] E E = C^e 7 , le a ,epj E B = le a ,epj E B =0, p E {e a )=p l a — : p E {e a )= — , 

and 

X E {x\y a )=y a e a , tt E {x\ y a ) = e a A e a + \cP ap y^ a A e~^, (4) 
(For more details, see [TTlliM] ). 

II. 3. Lie aflfgebroids. Let ta : A — > M be an affine bundle with associated vector bundle 
Ty : V — ► A/. Denote by t a+ : A + = Aff(A, K) ->• M the dual bundle whose fibre over x £ M 
consists of affine functions on the fibre A x . Note that this bundle has a distinguished section 
1_A £ F(t_4+ ) corresponding to the constant function 1 on A. We also consider the bidual bundle 
tj : A — > M whose fibre at x £ M is the vector space A x = (A x )* . Then, A may be identified 
with an affine subbundle of A via the inclusion ij± : A — > A given by i^(a)((p) — </?( a )i which 
is an injective affine map whose associated linear map is denoted by iy : V — » A. Thus, V 
may be identified with a vector subbundle of A. Using these facts, one can prove that there is 
a one-to-one correspondence between affine functions on A and linear functions on A. On the 
other hand, there is an obvious one-to-one correspondence between affine functions on A and 
sections of A + . 

A Lie affgebroid structure on A consists of a Lie algebra structure [•, -]y on the space r(ry) of 
the sections of Ty : V — > M, a K-linear action D : T(t_a) x r(ry ) — * r(ry) of the sections of A 
on r(ry) and an affine map pa ■ A — > TM, the anchor map, satisfying the following conditions: 

• D X \Y,Z} V = \D X Y,Z] V + \Y,D X Z] V , 

• D x+ yZ = D x Z+\Y,Z] v , 
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. D x (fY) = fD x Y + p A (X)(f)Y, 
for X G T(t a ), Y,Z G T(ry) and / G C°°(M) (see [12 [26]). 

If ([•) 'lv,D, pa) is a Lie affgebroid structure on an affine bundle A then (V, [-, -]y , py) is a Lie 
algebroid, where py : V — > TAf is the vector bundle map associated with the affine morphism 
p A : A -> TM. 

A Lie affgebroid structure on an affine bundle T4 : „4 — > M induces a Lie algebroid structure 
([') 'Ja' P l) on t ne bidual bundle A such that 1a G L(r_^+) is a f-cocycle in the corresponding 
Lie algebroid cohomology, that is, d A lA = 0. Indeed, if Xq G L(t^) then for every section X 
of A there exists a function / G C°°(M) and a section X G T(ry) such that X = /_X"o and 

PjifXo+X) = /aa(Xo)+M*)> 
[/Xo + X^Xo + r]^ = (py(X)( 5 )-py (?)(/) +fp A (X ){g) 

-gp A (X )(f))X + [X,Y] V + fD Xa Y - gD X(! X. 

Conversely, let (U,{-, -Ju , Pu) be a Lie algebroid over M and cf> : U — > R be a 1-cocycle of 
(t/, [•, such that </>|;y x ^ 0, for all x £ M. Then, A — </3 _1 {l} is an affine bundle over 

M which admits a Lie affgebroid structure in such a way that (U, [•, •]{/, pu) may be identified 
with the bidual Lie algebroid (A, [•, -\ Al p A ) to A and, under this identification, the 1-cocycle 
1.4 : A — > K is just <j>. The affine bundle r A '■ A — > M is modelled on the vector bundle 
ry : V — _1 {O} — > M. In fact, if «y : V — > {/ and i.4 : *4 — > C/ are the canonical inclusions, 
then 

i v o\X,Y} v = [i v oX,i v oYlu, i v oDxY = [i A oX,iv°Y\u, 
Pa(X) = pu{ia°X), 

for X, Y G L(ry) and X 6 T(t a ). 

Let T4 : .A — > M be a Lie affgebroid modelled on the Lie algebroid ry : 7->M, Suppose 
that (x 1 ) are local coordinates on an open subset U of M and that {eo, e a } is a local basis of 
sections of tj : A — > M in [7 which is adapted to the 1-cocycle I.4, i.e., such that l A (eo) = 1 
and l A (e a ) — 0, for all a. Note that if {e°, e Q } is the dual basis of {eo, e a } then e° = 1.4. Denote 
by (x 1 , y°,y a ) the corresponding local coordinates on A. Then, the local equation defining the 
affine subbundle A (respectively, the vector subbundle V) of A is y° = 1 (respectively, y° = 0). 
Thus, (x l 7 y a ) may be considered as local coordinates on A and V. 

The standard example of a Lie affgebroid may be constructed as follows. Let r : M — > R be a 
fibration and t^q : JV — » M be the 1-jet bundle of local sections of r : M — > R. It is well known 
that 71,0 : J 11 ?" — > M is an affine bundle modelled on the vector bundle tt = (%m)\Vt '■ V T ~ * M 
where W is the vertical bundle of r : M — > M. Moreover, if £ is the usual coordinate on R 
and 77 is the closed 1-form on M given by r\ — r* (<i£) then we have the following identification 
J l T = {v E TM/r](v) = 1} (see, for instance, OS ). Note that Vr = {v G TM/<q{v) = 0}. Thus, 
the bidual bundle J x r to the affine bundle ti.o : JV — > M may be identified with the tangent 
bundle TM to M and, under this identification, the Lie algebroid structure on ttm ■ TM — ► M 
is the standard Lie algebroid structure and the 1-cocycle lji T on ttm ■ TM — > M is just the 
1-form 17. 
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II. 4. The Lagrangian formalism on Lie affgebroids. In this section, we will develop a 
geometric framework, which allows to write the Euler-Lagrange equations associated with a 
Lagrangian function L on a Lie affgebroid A in an intrinsic way (see [26]). 
Suppose that (t^ : A —* M, ry : V — > M, ([•, -]y, D, pa)) is a Lie affgebroid over M. Then, the 
bidual bundle r j : A — ► M to A admits a Lie algebroid structure ([■, •] j,Pj) in such a way 
that the section I a of the dual bundle A + is a 1-cocycle. 

Now, we consider the Lie algebroid prolongation (T A A, [•, •]'-*, p T -~) of the Lie algebroid 
(A, [•, p^) over the fibration rx : A — ► M with vector bundle projection r- 4 : T A A — > *4. 
If (x l ) are local coordinates on an open subset U oi M and {eo, e Q } is a local basis of sections 
of the vector bundle ~~ * U adapted to 1a, then {To, T a , V a } is a local basis of sections 
of the vector bundle {r T £)- 1 {t^ 1 {U)) -> t^(U), where 

f (a) = (e Q (T A ( & )) lP ^^, T Q (a) = (^(^(aJJ.p^J, K(a) = (0, ^ ), (6) 

(x l , y a ) are the local coordinates on „4 induced by the local coordinates [x l ) and the basis {e a } 
and pg, are the components of the anchor map p j. We also have that 

[f , f a f£ = CJ Q T 7 , [T Q , fpf£ = C^f 7 , 

[To, = [T Q , = [V a , ^ - 0, (?) 

p r £(f ) = Pl-^r, P T £(f a ) = pi-^r, p r £(y a ) = JL, 

where Cq 1 ^ and Cf^g are the structure functions of the Lie bracket [•, with respect to the 
basis {e , e Q }. Note that, if {T °, f Q , K Q } is the dual basis of {f , f Q , V Q }, then T° is globally 
defined and it is a 1-cocycle. We will denote by 0o the 1-cocycle T°. Thus, we have that 

0o(a)(b,X a ) = 1.4(b), for (b,X a ) £ T/A 

One may also consider the vertical endomorphism S : T A A — > as a section of the vector 

bundle T-^^4 ® — > .A, whose local expression is (see [2"o] 1 

S =(T a -y a <t >0 )®V a - (8) 
A section £ of t- 4 : — > ^1. is said to be a second order differential equation (SODE) on A 

if MO = 1 and S £, = 0. If f € T^* 4 ) is a S0DE then C = To + y°T Q + £ a V a , where £ Q are 
local functions on A, and 

Now, a curve 7 : / C R ^ .4 in .4 is said to be admissible if 0^07 = {ja o 7) or, 
equivalently, (ij\(-f(t)),j(t)) £ T A t ^A, for all £ £ I, ij\ : A —> A being the canonical inclusion. 
We will denote by Adm(A) the space of admissible curves on A. Thus, if 7(4) = (x l (t) , y a (t)) , 
for all t £ I, then 7 is an admissible curve if and only if 

dx l 

-jjT =P l o + P l a V a , for *G{l,...,m}. 
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It is clear that if £ is a SODE then the integral curves of the vector field p T j^ (£) are admissible. 
On the other hand, let L : A — > K be a Lagrangian function. Then, we introduce the Poincare- 
Cartan l-section 0/, € T 
ciated with L defined by 



Cartan l-section O^ € r((T~ 4 )*) and £/ie Poincare-Cartan 2-section Ql E r(A (t- 4 )*) asso- 



6 L = L0 O + (d r ■*&) oS, fi L = -d r ^Gl. (9) 
From 0, © and ©, we obtain that 



n t = 



(ho(^(^)) - + cg a /) - a (io) 



d 2 L p 1 / i 8 2 L _ i d 2 L dL \ 

dy a dyP * + 2 V 13 dx l dy a Pa dx l dyP + dy^ af3 J 



> a A i 



where a = T a - y a $ ih ^ a = V a - Co and Co = T + y a T a + £g% is an arbitrary SODE. 
Now, a curve 7 : 7 = (— e, e) C R — > „4 in A is a solution of the Euler- Lagrange equations 
associated with L if and only if 7 is admissible and i(i A (-y(t)) .^(t))^L{l{t)) = 0, for all t. 
If 7(t) = (x 1 {t) , y a (t)) then 7 is a solution of the Euler-Lagrange equations if and only if 

for i £ {1, . . . , m} and a G {1, . . . , n}. 

If we denote by C(A + ) the set of curves in A + , we can define the Euler-Lagrange operator 
SL : Adm(A) -> C(-4+) by 

5i 7(t) (a) =n i (7(t))((u(7(t)),7(*)) > (S,« 7 (t))) s (12) 

for 7 € Adm(^4) and a, £ «4r^(7(t))i where v 7 ( t ) £ T^^A is such that (a, u 7 (t)) £ T^ t ^A, for all 
<. It is easy to prove that the map SL doesn't depend on the chosen tangent vector v^m. From 
(fTUl) . we deduce that its local expression is 

d 



where {e°, e a } is the dual basis of {eo, e a }. Then the Euler-Lagrange differential equations read 
as SL = 0. 

/ d 2 L \ 

The Lagrangian L is regular if and only if the matrix (W a g) = I — — —5- ) is regular or, in a 

\dy a oyP J 

intrinsic way, if the pair (f2z,, 4>o) is a cosymplectic structure on T A A, that is, 

{ fl L A . . . A»L A0 o }(a) ^0, d r ' A - A n jC , = and d^^o = 0, for all a e A 

(n 

Note that the first condition is equivalent to the fact that the map \>l '■ T A A — * (T A A)* defined 

by 

is an isomorphism of vector bundles. 
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Now, suppose that the Lagrangian L is regular and let Rl S r(r~ 4 ) be the Reeb section of the 
cosymplectic structure (TIlt&o) characterized by the following conditions 

iR L tt L = and i Rll <h = 1. 

Then, Rl is the unique Lagrangian SODE associated with L, that is, the integral curves of the 
vector field p t -~(Rl) are solutions of the Euler-Lagrange equations associated with L. In such 
a case, Rl is called the Euler-Lagrange section associated with L and its local expression is 

/ fit 8 2 t Fit \ 

R L =f + y«f a + ^[p^— (pi + V^tow + (Cfc + V*Cfo) w )v a , 
where (W al3 ) is the inverse matrix of (W a p). 

II. 5. The Hamiltonian formalism. In this section, we will develop a geometric framework, 
which allows to write the Hamilton equations associated with a Hamiltonian section on a Lie 
affgebroid (see [HCI5]). 

Suppose that (t.4 : A — ► M, ry : V — > M, ([•, D, p_^j) is a Lie affgebroid. Now, consider the 
prolongation (T A V* , [•, 'YjiPj) °f the bidual Lie algebroid (A, [-, -]^, pjj) over the fibration 
t^:V* -> M. 

Let be local coordinates on an open subset U of M and {eo, e^} be a local basis of sections 
of the vector bundle (U) — > {7 adapted to 1^ and 

[e ,ej_4 = C y Q e 7 , [e„, e^lj = C^e^ pj{e ) = Pl^-, P^( e ")=Pa^— ■ 

Denote by (x l ,y°, y a ) the induced local coordinates on A and by {x l ,yo, y a ) the dual coordinates 
on the dual vector bundle T4+ : A + — > M to „4. Then, (x l ,y a ) are local coordinates on y* and 
{eo, e Q , e Q } is a local basis of sections of the vector bundle r-' : T A V* — > V* , where 



.4 

e ' ' ' ' « - / 



>(*) = (-oWWU^J, = (^(^^ J, = (o, JL J. 



Using this local basis one may introduce local coordinates (x % , y a ; z°, z a ,v a ) on T A V* . A direct 
computation proves that 

[ g o, epf- = C y /9 e 7 , [e«, e fj = C2 /3 e 7 , 

[eo.ejj = [e a) e^]J = [e^e^J = 0, 

for all a and /3. Thus, if {e°, e Q , e Q } is the dual basis of {eo, e Q , e a } then 

= -\cr Qa e° Ae<* -\clfe a Ae", d r ^* e° = V = 0, 

for / e C°°(V*). 
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Let fi : A + — > V* be the canonical projection given by fi((p) — cp l , for G .A+, with a; G M, 
where ip l G V^* is the linear map associated with the affine map ip and h : V* —> A + be a 
Hamiltonian section of /!. 

Now, we consider the Lie algebroid prolongation T A A + of the Lie algebroid A over T4+ : 
.4+ — > M with vector bundle projection Tj" 4+ : T- 4 4 + — > .4+ (see Section fll. 2|) . Then, we 

may introduce the map Th : T A V* -> defined by Th(&,X a ) = (&,(T a h)(X a )), for 

(a, X Q ) G T A V*, with a G V*. It is easy to prove that the pair (Th,h) is a Lie algebroid 

morphism between the Lie algcbroids t t -? : T A V* — > V* and r~- 4+ : T- 4 ^ — » _4 + . 

.A ,A 

Next, denote by A h and the sections of the vector bundles (T A V*)*~ -> V* and A 2 (T A V*)* -> 
V* given by 

A h = (T/i,fc)*(Aj) ) O ft = (T/ l) /i)*(^), (13) 
where A_^ and Ojj are the Liouville section and the canonical symplectic section, respectively, 
associated with the Lie algebroid A. Note that Clh = —d rAy Xh- 
On the other hand, let 7? : T A V* — > R be the section of (T A V*)* -> ^* defined by 

7? (a,X I/ ) = U(a), (14) 

for (a, X v ) G T A V*, with veF. Note that if pn : T- 4 ^* — > A is the canonical projection on 
the first factor then (pri,T v ) is a morphism between the Lie algebroids r-' : T A V* — > V* and 
rr : -4 — > M and (p^i,Ty)*(l^i) = 77. Thus, since 1^4 is a 1-cocycle of t? : A — > M, we deduce 
that 77 is a 1-cocycle of the Lie algebroid t t £ : T A V* — > V*. 

Suppose that h{x l ,y a ) = {x 1 , —H(x^,yj3),y a ) and that {e°,e Q ,e Q } is the dual basis of {eo,e a , 
e a }. Then 77 = e° and, from (UJ, ([13]) and the definition of the map Th, it follows that 

1 /9 TJ flTJ 

n h = e a A e<* + -C2 y 7 e a A e' 3 + (pj,^ - Cg^S" A e° + — e Q A g°. (15) 

Thus, it is easy to prove that the pair (0,^,7]) is a cosymplectic structure on the Lie algebroid 

T r y ■ t a v* — > 

Remark II. 1. Let T^V* be the prolongation of the Lie algebroid V over the projection 
Ty : V* — > M. Denote by Ay and Sly the Liouville section and the canonical symplectic 
section, respectively, of V and by (iy,ld) '■ T V V* — » T-^V* the canonical inclusion. Then, 
using ^ , <P~S]) , (TT4^ and the fact that /i o h = Id, we obtain that 

(iv,Id)*(\ h ) = \ v , (i v ,ld)*(r])=0. 

Thus, since (iy, Id) is a Lie algebroid morphism, we also deduce that 

{i v ,id)*(n h ) = Q v . 



Now, let Rh G Y{t~) be the Reeb section of the cosymplectic structure (f^,??) characterized 
by the following conditions 

iR h ^h = and i Rh i] = 1. 
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With respect to the basis {eo, e a , e, 



Rh = e + q — 




(16) 



Thus, the vector field p T j (Rh) is locally given by 



- « + + ( - *s? + + ^ssO %: (17) 



and the integral curves of Rh (i.e., the integral curves of p-' (Rh)) are just the solutions of the 
Hamilton equations for h, 




for i £ {1, . . . , m} and a G {1, . . . , n}. 

Next, we will present an alternative approach in order to obtain the Hamilton equations. For 
this purpose, we will use the notion of an aff-Poisson structure on an AV-bundle which was 
introduced in [5] (see also [TO]). 

Let tz : Z — > M be an affine bundle of rank 1 modelled on the trivial vector bundle tmxr : 
AI x E — > M, that is, r z : Z — > M is an AV-bundle in the terminology of [TO] , 
Then, we have an action of K on the fibres of Z. This action induces a vector field JT^ on Z 
which is vertical with respect to the projection tz : Z — > A/. 

On the other hand, there exists a one-to-one correspondence between the space of sections of 
tz : Z — > M, r(rz), and the set 



In fact, if /i G rfrz) and (x l , s) are local fibred coordinates on Z such that Xz — tt and ft is 
locally defined by h(x % ) — (x z , —H(x 1 )), then the function Fh on Z is locally given by 



{F h G C°°(Z)/X Z (F A ) 



1}. 



d 




(18) 



{ v }:r(r z )xr(r^rM 



which satisfies the following properties: 

i) Skew-symmetric: {hi, h%] = — {h>2, hi}. 

ii) Jacobi identity: 



{hi, {h 2 , h 3 }} v + {h 2 , {h 3 , hi}} v + {h 3 , {hi, h 2 }}v = 0, 



where {■, -}y is the affine-linear part of the bi-affine bracket, 
iii) If h G r(rz) then 



{h,-}:T(T Z )^C°°(M), h'^{h,h'}, 



is an affine derivation. 
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Condition iii) implies that, for each h £ T(t z ) the linear part {h,-} v ■ C°°(M) -> C°°(M) 
of the affine map {h,-} : T(tz) — * C°°(M) defines a vector field on AI, which is called the 
Hamiltonian vector field of h (see |10j). 

In [TU], the authors proved that there is a one-to-one correspondence between aff-Poisson bra- 
ckets {•, ■} on tz : Z — > M and Poisson brackets {•, -}n on Z which are A^-invariant, i.e., which 
are associated with Poisson 2-vectors II on Z such that £jc z II = 0. This correspondence is 
determined by 

{hi,h 2 }oT Z = {F hl ,F h2 }n, for h 1 ,h 2 £T(T Z ). 
Note that the function {Fh ± , F/ l2 }u on Z is r^-projectable, i.e., Cx z {Fhi > Fh 2 }n = (because 
the Poisson 2-vector IT is X^-invariant). 
Using this correspondence we will prove the following result. 

Theorem II. 2. [12) Let T4 : A — > M be a Lie affgebroid modelled on the vector bundle 
Ty : V — > M. Denote by : A + — » M (resp., t v : V* — > M) the dual vector bundle to A 
(resp., to V) and by fi : A + — > V* the canonical projection. Then: 

i) fi : A + —> V* is an AV-bundle which admits an aff-Poisson structure. 

ii) If h : V* — > A + is a Hamiltonian section (that is, h £ then the Hamiltonian 
vector field of h with respect to the aff-Poisson structure is a vector field on V* whose 
integral curves are just the solutions of the Hamilton equations for h. 

Proof, i) It is clear that fx : A + — > V* is an AV-bundle. In fact, if a + £ A%, with x £ M, and 
t £ R then 

a+ + t = a+ + tl A (x). 

Thus, the /i-vertical vector field on A + is just the vertical lift 1^ of the section 1^ £ 

T(t_4+). Moreover, one may consider the Lie algebroid tt : A = (^4 + )* — > M and the corres- 
ponding linear Poisson 2-vector II4+ on A + . Then, using the fact that I.4 is a 1-cocycle of 
rr : A = (A + )* — > A/, it follows that the Poisson 2-vector n_4+ is X4+ -invariant. Therefore, 
IT_4+ induces an aff-Poisson structure {•, •} on // : .4+ — > U* which is characterized by the 
condition 

{hi,h 2 }o /j, = {F hl ,F h2 }n A+ , for h u h 2 £T(jx). (19) 
One may also prove this first part of the theorem using the relation between special Lie affgebroid 
structures on an affine bundle A' and aff-Poisson structures on the AV-bundle AV((A')$) (see 
Theorem 23 in [TU]). 

ii) From (JTHJ and JU), we deduce that the linear map {h, -}v ■ C°°(V*) -> C°°(y*) associated 
with the affine map {h, ■} : T(fi) — > (that is, the Hamiltonian vector field of ft.) is given 

by 

{/v}vM°M={^°mK + , for ^eC°°(y*). (20) 
Now, suppose that the local expression of h is 

h(x i ,y a ) = (x i ,-H{x s ,y p ),y a ). (21) 
On the other hand, using (J5]), we have that 

„ ,-9 9 ,9 9 _~ 9 9 1- 9 9 

^ + = *W A 9^ + Aft? A W a ~ C °^Wo A W a ~ * ft* (22) 
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Thus, from fZU\i , (|2"T|) and (j2"2"|) , we conclude that the Hamiltonian vector field of h is locally 
given by 

, i ,■ dH . d j , dH ^ dH \ d 

which proves our result (see (fl7|)). □ 

II. 6. The Legendre transformation and the equivalence between the Lagrangian and 
Hamiltonian formalisms. Let L : A — > R be a Lagrangian function and Ol G r((r~ 4 )*) be 
the Poincare-Cartan 1-section associated with L. We introduce the extended Legendre transfor- 
mation associated with L as the smooth map LegL : A — > A + defined by Legi,(&)(b) = 0i(a)(z), 
for a, b £ A x , where z £ T A A is such that pr\{z) — i^.(b), pri : T A A — > .A being the restriction 
to T-^A of the first canonical projection pr\ : A x T*4 — > A The map Legz, is well-defined and 
its local expression in fibred coordinates on A and A + is 

Le 9L (x\y") = ( X \L-—y«,^). (23) 

Thus, we can define the Legendre transformation associated with L, legL ■ A — » V*, by Zeg^ = 
/i o LegL- From (f23|) and since /z(a; 1 , yo, y a ) = {x l , we have that 

leg L (x\y a ) = (x\ — ). (24) 

The maps Leg_L and /e#L induce the maps TLegL ■ T A A — > T A A + and TlegL ■ T A A — > T- 4 ^* 
defined by 

(TLe.g L )(b,X a ) = (b, (T a Le 5i )(X a )), (TZe ffL )(b, X a ) = (b, (T a 2e 5L )(X a )), (25) 
for a £ A and (b, X a ) e T/A 

Now, let {Ta,T a ,V a } (respectively, {eo, e a , e~o, e a }) be a local basis of T(t~^) as in Section 
QT4] (respectively, of T(t^ 4+ ) as in Section Hi. 2p and denote by (x l , y a ; z°, z a , v a ) (respectively, 
(x l , yo, y a 'i z°, z a , vo, v a j) the corresponding local coordinates on T A A (respectively, T A A + ). 
In addition, suppose that (x 1 , y a ', z°, z a , v a ) are local coordinates on T A V* as in Section ["11.51 
Then, from (f2"3"|) . (12"4")) and ([23]) . we deduce that the local expression of the maps TLegL and 
TlegL is 

TLeg L (x\ y a ;z°, z a ,v a )=(x\L - ^y a , ^;z°, z a ,z°pl(^ - ^^y a ) 

+Z P «W dx*dyP V } V dy°dy0 y ' Z Po dx*dy" ( ' 

+Z P0 dx*dy a +V dy^dyP' 1 



Tleg L (x\y«;z°,z«,v*) = (a* "- + ** p <, 



,/3 



' ' ^dx^y 01 rp dx^dy a 
d 2 L , 



(27) 



dy a dyV 

Thus, using ([I]), JTOJ), (ESJ) and (f2l)|) . we can prove the following result. 
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Theorem II. 3. |12j The pair (TLegL, LegL) is a morphism between the Lie algebroids (T A A, 
[•, -Yj~ , pj) and (T A A + , p T j + ). Moreover, if Ql and fir, (respectively, Xj and Q^) 
are the Poincare-Cartan 1-section and 2-section associated with L (respectively, the Liouville 
1-section and the canonical symplectic section associated with A) then 

(TLeg L ,Leg L )*{\ I ) = Q Ll (TLeg L , Leg L )* (%) = fl L . (28) 
From (ED, it follows 



Proposition II. 4. |12j The Lagrangian L is regular if and only if the Legendre transformation 
legL '■ A — » V* is a local diffeomorphism. 

Next, we will assume that L is hyperregular, that is, legL is a global diffeomorphism. Then, from 
(|25p and Theorem III. 3( we conclude that the pair (TlegL, legL) is a Lie algebroid isomorphism. 
Moreover, we can consider the Hamiltonian section h : V* —* A + defined by 

h = Leg L o leg^ 1 , (29) 

the corresponding cosymplectic structure (Qh,rf) on the Lie algebroid t~ : T A V* — > V* and 
the Hamiltonian section Rh S T(t t ~). 

Using ([15 ]) . (|37 ]) . ([25 ]) . ([21?] ) and Theorem |IL3] we deduce that 

Theorem II. 5. |12j // the Lagrangian L is hyperregular then the Euler- Lagrange section Rl 
associated with L and the Hamiltonian section Rh associated with h satisfy the following relation 

Rh ° legL = TlegL ° Rl- 

Moreover, if 7 : / — > A is a solution of the Euler-Lagrange equations associated with L, then 
legL 7 ■ I - * V* is a solution of the Hamilton equations associated with h and, conversely, if 
7 : / — > V* is a solution of the Hamilton equations for h then 7 = leg^ 1 07 is a solution of the 
Euler-Lagrange equations for L. 

Note that if the local expression of the inverse of the Legendre transformation legL ■ A — > V* 
is 

leg^ix^yp) = (x\ y a (x j , y )), 
then the Hamiltonian section h is locally given by 

h(x\y a ) = (x\-H(x j ,yp),y a ), 

where the function H is 

H{x\y a ) = y f} y (x i ,y a )- L{x\y {x\y a )). 

Thus, 

— =-— — = (30) 

dx l \(xi, W ) dx l \ (xi ,y<*{xi ^p)) 9 V0\( x i,y a ) 



Therefore, 



legl 1 {x i ,y a ) = (x\—). (31) 
oy a 
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III. AFFINELY CONSTRAINED LAGRANGIAN SYSTEMS 

We start with a free Lagrangian system on a Lie affgebroid A of rank n. Now, we plug in 
some nonholonomic affine constraints described by an affine subbbundle B of rank n — r of the 
bundle A of admissible directions, that is, we have an affine bundle rg : B — > M with associated 
vector bundle tjj b : Ub — > M and the corresponding inclusions ig : B A and ijj B '■ V, 
Ty : V — > M being the vector bundle associated with the affine bundle T4 : A — > M. If we 
impose to the admissible solution curves 7(f) the condition to stay on the manifold B, we arrive 
at the equations 5L 7 ^ = X(t) and 7(f) G B, where the constraint force A(t) G ^S- A { (t)) ^ s ^° ^ e 
determined. In the standard case (.A = J 1 r, r : M — > R being a fibration), A takes values in the 
annihilator of the constraint submanifold B. Therefore, in the case of a general Lie affgebroid, 
the natural Lagrange- d'Alembert equations one should pose are 

8L l[t) G B° Bh(t)) and j(t) G B, 

where B° — {ip G A + /(pm = 0} is the affine annihilator of B which is a vector subbundle of A + 
with rank equal to r. 

In more explicit terms, we look for curves 7(f) € A such that 

- 7 G Adm(A), 

- 7(*) e B T^( 7 (t)), 

- there exists X(t) G ^°AMt)) snc ^ tnat <^7(t) = M*)- 

If 7(f) is one of such curves, then (7(f), 7(f)) is a curve in T A A with i(-y(t)Mt))^Q = 1- More- 
over, since 7(f) is in B, we have that 7(f) is tangent to B, that is, (7(f), 7(f)) G T A B with 
*(7(i),7(t))^o ~ 1- Note that, since B is an affine subbundle of A, the Lie algebroid prolongation 

TB . _^ g of algebroid ^ 

over the fibration Tg : B -> M is well-defined. Under 
some regularity conditions (to be made precise later on), we may assume that the above admi- 
ssible curves are integral curves of a section X which we will assume that it is a SODE section 
taking values in T A B. Then, from (fT2|) . we deduce that 

(*xfi£)(b) e K, for all b G B, 

where B° is the vector bundle over B whose fibre at the point b G B is 

B£ = {a e (T^Ay/a(b',Vb) = 0,V(b> b ) G T^A such that b' G B TA(h) }. 

Based on the previous arguments, we may reformulate geometrically our problem as the search 
for such a SODE X (defined at least on a neighborhood of B) satisfying ix(b)^i(b) G 
*x(b)^o(b) = 1 and X(b) G T^B, at every point b G B. 
Now, we will prove the following result. 

Proposition III.l. If S is the vertical endomorphism then the vector bundles overB, B° — > B 
and S*{T A B)° — > B, are equal. In other words, 

B,° = S*{T h J B)°, for all b G B. 

Proof. First, we will see that the map 

S* : (T h A B)° -> S*(T h X B)° 
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is a linear isomorphism. 

In fact, suppose that a e (T^B)° and S*a = 0. Then, we have that 

a(0, Vh) = 0, for all Vb & T^A such that (TbT^)^) = 0. 

Moreover, if (a,itb) £ ^b* 4 ^ then, since T^t^ '■ T^B — » T TA i^M is a linear epimorphism, it 
follows that there exists i>b S Tb£> satisfying 

(7bT4)(u b ) = ^(5). 

Thus, using that a G we deduce that 

a(a, Ub) = £*(5, «b) + ck(0, Ub — «b) = 0. 

Therefore, a = 0. 

This implies that S* : (T h A B)° — > S*(T h A B)° is a linear isomorphism and 

dimS*(T^B)° = r. 

On the other hand, since 7b T4 : TbS — ► T^^M is a linear epimorphism, we obtain that 

dimBl =r = dim S*(T h A B)° . 
Finally, we will prove that S*(T A B)° C B°. 

In fact, if a e (T b J B)°, (b', w b ) £ T h J A and b' £ B TB(b ), it follows that 

S(b',Vb) = (0,<), with<GT b B. 

Consequently, 

(S*a)(b',v h ) = 0. 

□ 

Proposition IIII. ll suggests us to introduce the following definition. 

Definition III. 2. A nonholonomically constrained Lagrangian system on a Lie affgebroid A 
is a pair (L,B), where L : A — ► R is a C°° -Lagrangian function, and i& : B ^ A is a 
smooth affine subbundle of A, called the constraint affine subbundle. By a solution of the 

nonholonomically constrained Lagrangian system (L,B) we mean a section X £ ^Xt- 4 ) which 
is a SODE on B and satisfies the Lagrange- d'Alembert equations 

{ix^L)\B G r(r st(T ^ B)0 ), 

{ix<j>a)\B = 1, (32) 

X\B g r(rj), 

where T s *^q-A^o is the vector bundle projection of S*(T A B)° — > B. 

With a slight abuse of language, we will interchangeably refer to a solution of the constrained 
Lagrangian system as a section or the collection of its corresponding integral curves. 
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Remark III. 3. We want to stress that a solution of the Lagrange-d'Alembert equations needs 
to be defined only over B, but for practical purposes we consider it extended to A (or just to 
a neighborhood of B in A). We will not make any notational distinction between a solution 
on B and any of its extensions. Solutions which coincide on B will be considered as equal. In 
accordance with this convention, by a SODE on B we mean a section of T A B which is the 
restriction to B of some SODE defined in a neighborhood of B. (} 

Now, we will analyze the form of the Lagrange-d'Alembert equations in local coordinates. Let 
(x z ) be local coordinates on an open subset U of M and {eo, e Q } be a local basis of sections of 
the vector bundle ^(U) — > U adapted to 1^. Denote by (x l ,y°,y a ) the corresponding local 
coordinates on A and suppose that {/J^e a } is a local basis of the vector bundle U% — > M and 
that e® is a local section of the affine bundle B — > M such that eo~ e® = age Q . Then, using that 
an element eo(x) + y a e a (x) of A x belongs to B x if and only if (eo(x) + y a e a (x)) — e® (x) G (Ub) x , 
we deduce that the local equations defining the constrained subbundle B as an affine subbundle 
of A are 

* a = + /0° = 0, fOT a11 a G {1, . . . , r}, 
where /ig = ^afi. Observe that d T ' A - 4 * a G r((rj*)*) verifies that 

(d T ^* a )(b)(a,Jf b ) = 0, for all (5, X b ) G t/b, 

and, therefore, (d T ' 4 - 4 * a )(b) G (T b J B)°, for every b G B and a G {1, . . . , r}. 
On the other hand, using (TT|) and (J7J), we deduce that 

{0o, r Q , being the dual basis of {f , f a , V a }. 

Thus, since the matrix has maximum rank r, it follows that the sections {d 7 " A " /l Vl /a }a=i....,r 
are linearly independent and, using that rank(T A B)° = r, we deduce that {(d T A ^ a )\ B }a=i,...,r 
is a local basis of sections of (T A B)° — > B. 

Therefore, using © and the fact that S* : (T^B)° -> S*(T^B)° is an isomorphism, for all 
b e 6, we obtain that {/x°6 Q! } a= i ! ... jr is a local basis of sections of 5 , *(T- A B)° — > B, where 
e a = f a - y a 4> . 

Consequently, the constrained motion equations (|32p can be written as 

ix4>o = 1, (33) 

along the points of B, where A a are some Lagrange multipliers to be determined. 
Now, a section X G T(r^ A ) is of the form X = z°f + z a f a + v a V a . If we assume that the 
Lagrangian L is regular, then the first two equations in (|33[) imply that any solution X has to 
be a SODE, that is, X = T + y a T a + X a V a and, then, the first equation in ([3"3"|) becomes (see 

0330 

f) 2 T r) 2 T BT BT 
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Thus, if Rl is the Euler-Lagrange section associated with L, X = Rl — W a ^X a fi'^V a and then, 
the third equation in (|33[) implies that 



(d 1 ^ a )(R L ) + h(d J A ^ a )(-W a "fi b V a ) = O, for all ae{l,...,r}. 

As a consequence, we get that there exists a unique solution of the Lagrange-d'Alembert equa- 
tions if and only if the matrix 



c ab (b) = (-w^<)(b) 



(34) 



is regular, for all b € 23. 

From (|3"3"1) , we deduce that the differential equations for the integral curves of the vector field 
p T ~(X) are the Lagrange-d'Alembert differential equations, which read 



d_ f dL 
~dt \dy* 



~dt 
dL 



p l + p l «y a , 



dL 



(35) 



Mg + M«J/ Q = 0, 
with i G {1, . . . , m}, a S {1, . . . , n} and a G {1, . . . , r}. 



Remark III. 4. In the above discussion, we obtained a complete set of differential equations 
that determine the dynamics. Now, we will analyze the form of the Lagrange-d'Alembert 
equations in terms of adapted coordinates to B. Consider local coordinates (x l ) on a open set 
U of M and the local basis {<j>A~}A=i,...,n-r °f sections of C/g. Complete it to a local basis of 
sections {eo, 4>a, 4>a} of the vector bundle T j^~{U) — > U adapted to 1a- Then, in coordinates 
(x z ,y°, y A ,y a ) adapted to this basis, the equations defining the constrained subbundle B as an 
affine subbundle of A (respectively, as an affine subbundle of A) are y a = (respectively, y° = 1 
and y a = 0). Thus, the Lagrange-d'Alembert differential equations read as 



dx l 

-dt =P ° 



PaV 



d dL , 
di^dy^' 



Pa 



dL 



(C 7 



AO 



CIbV 1 



dL 



0. 







IV. Solution of Lagrange-d'Alembert equations 

In this section, we will perform a precise global analysis of the existence and uniqueness of the 
solution of Lagrange-d'Alembert equations. In what follows, we will assume that the Lagrangian 
L is regular at least in a neighborhood of B. 

Definition IV. 1. A constrained Lagrangian system (L, B) is said to be regular if the Lagrange- 
dAlembert equations have a unique solution. 
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In order to characterize geometrically these nonholonomic systems which are regular, we define 
the vector subbundle F C T A A\ B -> B whose fibre at point b e B is F h = \>1 1 {S*{% A B)°), 
where \>l : T A A — > (T" 4 ^)* is the vector bundle isomorphism defined by 

b L (AT) = i x £l L + 0o(^)0o, for all X e T^A (36) 

More explicitly, 

Fb = {l£ A/ exists X a s.t. \> L (X) = Xa» a J a \b}- 
From the definition, it is clear that the rank of F is rank(F) = rank(T A B)° — rank(T A A) — 
rank(T A B) = rank(B°) = r. If we consider the sections Z a e r( T j 4 ) such that \>h{Z a ) = fJ>a& a S 
with a 6 {I, . . . , r}, then {Z } is a local basis of sections of F. Moreover, if Rl is the Euler- 
Lagrange section associated with L, we have that 

(i Za Q L )(R L ) + Mz a )MRL) = n a J a {R L ) = o 

which implies that cf>o(Z a ) = 0. Therefore, Z a is completely characterized by the conditions 

iz a M L = [i a J a and i Za <t>o = 0. 

In addition, using these two conditions, we conclude that the local expression of Z a is the 
following 

Z a = -W a ^%V a , (37) 
where (W a ^) is the inverse matrix of (W a p)- Thus, the matrix defined in (l34|) is just C ab = 
(d TAA V a ){Z b ) = -W al3 Li a 0f i h a and we will denote it by C = (C ab ). 

A second important geometric object is the vector subbundle G C T A A\& — » B whose annihi- 
lator at a point b S B is S*{T^B)°. We can consider the subbundle G 1 - C F A A\ B — > B, the 
orthogonal to G with respect to the cosymplectic structure (f^z,, 0o), which is given by 

G[ = {le T^A/(i x Si L + <h(X)<h)\G b = 0}, for b e B. 
Note that G^ = t>7/(G£) = F b , for all b e B. Moreover, we obtain 

Proposition IV. 2. Gb is coisotropic in (T A A, fii,(b), 0o(b)), /or a// beS. in otter words, 
Gi C G b , /or b e B. 

Proof. In fact, using ([3T|) and since i* 1 (respectively, G°) is locally generated by {2' a }a=i,...,r 
(respectively, {5*(d T ^' 4 * a )} a= i ! ... ir ), we deduce that 

Gi =F h C G b , for all b 6 B. 

□ 

Now, we introduce the section A^ of the vector bundle A 2 (T A A) — > .4 defined by 

A i (a,/3) = L (b^ 1 (a)>i 1 (^)) ! (38) 

for a,/J £ (T" 4 ^)*. Ai is the (algebraic) Poisson structure (on the vector bundle T A A — » .4) 
associated with the cosymplectic structure (Ql,4>o)- We will denote by $\ L : (7^.4)* — > T A A 
the vector bundle morphism given by 

U L (a) = i a A L , iorae(T A A)*. 
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Note that 

MfoA<*)) =n L Q>l 1 (a),R L ) =0. 
On the other hand, it is clear that 

a = hi^ia)) = h-\ a) ^L + ct>ofyl,H a ))4>o 

and, thus 

a(R L ) = MK\a))- (39) 

Moreover, from ((551) , we have that 
which implies that (see ([39])) 

P{Ut{ot)) = P{-\>l 1 {a) + a{R L )R L ). 

Therefore, we have proved that 

(t Ai (o) = -\>l\a) +a(R L )R L , for a e {T A A)*. (40) 

Next, we consider the vector subbundle G ±,Al - C T A A\& — > £>, the orthogonal to G with respect 
to the Poisson structure (A^)^, which is given by 

Gi AL =h L (G° h ), forbeS. 

From (gO]) and since S*(d TAA ^ a ){R L ) = 0, we deduce that 

U L (s*{d TAA -* a )\ B ) = -z a . 

Consequently, we obtain that 

Gb Al = Gb = F b , for b e B. 

Now, we will consider the vector subbundle H C T A A\& — » £> whose fibre at point b 6 B is 
i?b = T^B (~1 Gb- We will denote, as above, by H ± the orthogonal to H with respect to the 
cosymplectic structure (fix,, 0o) and by H ± ' Al the orthogonal to -ff with respect to the Poisson 
structure A^. 

Using ([8]) and (fTO]) . we have that is&L = 0, that is, 

Cl L (SX,Y) + Q L (X,SY) = 0, foraU ijerf^), 

or, equivalently, 

isxto L = -S*(i x Sl L ), for all 16^). (41) 

Next, denote by grad 1 $> a the gradient section corresponding to the function ^f a with respect to 
the cosymplectic structure (0/,, cf>o), that is, grad = b^ 1 {d TAA ^ a ). Then, using (j4*Tj) and the 
fact that S*4>o = 0, we have that i^sgrad *°)| B ^ i = ~*.z a ^L an d 4>o{S 9 ra d ^ a ) = 4>o{Za) = 
which implies that 

(5 S rad * Q )| e = -Z„, for all a e {!,..., r}. (42) 
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Since the sections {d rJ l - 4 \I> a } are linearly independent, it follows that the sections {grad ^ a } 
also are linearly independent. Thus, using (|4*2"|) . we conclude that the sections {(grad x E' a )|£, Z a } 
are linearly independent. On the other hand, 

H° = (T A B)° + G° (43) 

and, therefore, the space of sections of H 1 ^ is locally generated by {(grad ^ a )|g,Z a }. This 
implies that {(grad v E ,a )|g, Z a } is a local basis of sections of H and rank(H ) = corank(H) = 
2r. 

Now, using we deduce that 

h ±,Al = ( T AtyX,A L + p 

On the other hand, we introduce the Hamiltonian section X^i associated with the function *S> a 
with respect to which is defined by 

= ~h L (d T * A * a ), for a G {1, ... , r}. (44) 

From (|4U|) . we have that 

h(X&) = d T * A * a - (d TXA * a )(R L )4> . (45) 

Thus, using that the sections {d TAA ^! a , 4>q} are independent, it follows that the sections {X^i} 
are also independent. Moreover, from (|4"5|) . we obtain that 

X$i = grad f a - (d 7 * A * a )(R L )R L 

which implies that 

SX*i = -Z a , fcroe{l r}. 

Therefore, {X^i , Z a } is a local basis of sections of H ±,Al and rank(H A -' AL ) = corank(H) = 2r. 
The relation among these objects is described by the next result. 

Theorem IV. 3. The following properties are equivalent: 

(1) The constrained Lagrangian system (L, B) is regular, 

(2) the matrices C = (C ab ) are non- singular, 

(3) T A B n F h = {0}, for all b G B, 

(4) H h n H£ = {0}, for all beB, 

(5) H h n H^ Al = {0}, for all beB. 

Proof. [(1) 44> (2)] This result was proved in Section HTT1 

[(2) <=> (3)] Suppose that C is non-singular and let be X G 7"/,B n F b . Thus, X = 

Ea=i A ^a(b) and (rf^tf"^)^) = 0, for all a G {1, ... , r}, which implies that 

J] \ c (d TAA * a )(b)(Z c (b)) = £ A c C ac = 0. 

c=l c=l 

Therefore, we deduce that A c = 0, for all c, and consequently X = 0. 

(<=) Conversely, take an arbitrary linear combination of columns of C at some point b such that 

r 

A c C ac (b) = 0, for all a G {1, ... , r}. 

c=l 
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Thus, X)c=i ^cZ c (b) E T^B which implies that Y^c=i ^cZ c (b) = 0, and hence A c = 0, for all 
ce{l,...,r}. 

[(3] <=> (4)] (=^) Let b E B, then we have that T^Bn F b = T^B n G£ = {0} and T h A A = 
T^B 8 Gt (note that dimT^B + dimG^ = dimT^A). Hence, from Proposition [TvT2l we 
obtain that 

G b = {T h J B n G b ) © Gi = H h Q) Gb • 

If X G i/b H i?^ then, from the above decomposition, we also have that X G G^ and, thus, 
X = 0. 

O) Conversely, let b G 6 and take an element X G 7^6 nFb = T^BdG^ C 7^6 n G b = # b . 
Since, ft L (b)(X,Y~) = and <p (b)(X) = 0, for all Y E iJ b , we conclude that X E and, 
therefore, X = 0. 

[(4) <^> (5)] (=>) Let b G B, then we have that H\> H — {0} and, as a consequence, the 
restriction (f^^ b , 4>Q h ) to -ffb of the cosymplectic structure (fi^, </>o) is a cosymplectic structure 
on i?b- In other words, the map 

bf b : H h -> tf b *, X G Jf b ^ (» x O L (b) + o (b)(X)<Mb)V b e #b, 

is a linear isomorphism. Thus, we can consider the Reeb vector R H (b) of the cosymplectic 
vector space (-ffb, ^£ b 5 ^o^) which is characterized by the following conditions 

i R H {h) n^ = and i RH(h) 4>^ = 1. (46) 

Now, we will prove that H h n %a l = {0}. 

Suppose that a G -ff b anc ^ ttA L (a) G -ff b . Then, using (JIU]) and the fact that </>o(b)(jiA i (a)) = 0, 
we deduce that 

i w{a) Sl L {b) = -a + a(R L (b))M^)- (47) 
Therefore, from (|46|) and since a G i7 b an d tlAi,( a ) G -ff b , we have that 

Q = (i tALia) n L (b))(R H (b))=a(R L (b)) 
which implies that (see (14T10 

^ b (tt Ai («)) = o. 

Consequently, JJa^ (a) = 0. 

(4=) Let b G B and take an element X E H h n H^. Then 

(ixfii(b) + <Mb)(X)0 o (b))| Hb = 0. 
Thus, using that X E i? b , it follows that 

<p (b)(X) = and a = i x tt L (b) E H£. 
Therefore, a(R L (b)) = and a = \>l(X). This implies that (see (JM])) 

X = -j| At (a)GtlA t ( J ff b )=if b L ' Ai . 

Consequently, X = 0. 

□ 

Proposition IV. 4. Conditions (3), (4) and (5) in Theorem \ I V.3\ are equivalent, respectively, 
to 
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(3') T A A lB = T A B®F, 
(4 ! ) T-*A\b = H®H ± , 
(5') T A A\ B = H ®H J - Al -. 

Proof. The equivalence of (3) and (3') follows by computing the dimension of the corresponding 
spaces. The ranks of T A A, T A B and F are 

rank(T A A) = 2 rank(A) + 1, 
rank(T A B) = rank(A) + rank(B) + 1, 
rank(F) = rank(B°) = rank(A) — rank(B). 

Thus, rank(T A A) = rank(T A B) + rank(F), and the result follows. The equivalence between 
(4) and (4') is obvious, since we are assuming that the free Lagrangian is regular, i.e., (fix,, <fio) 
is a cosymplectic structure on T A A. Finally, the equivalence of (5) and (5') is also obvious 
because we have that rank(H^-> AL ) = corank(H). 

□ 

IV. 1. Projectors. Now, we can express the constrained dynamical section in terms of the free 
dynamical section by projecting to the adequate space, either T A B or H, according to each of 
the above decompositions of T A A\$. Of course, both procedures give the same result. 

IV. 1.1. Projection to T A B. Assuming that the constrained Lagrangian system is regular, we 
have a direct sum decomposition 

T h A A = T h A B®F h , for all beB, 

where we recall that the subbundle F C T A A is defined by F = b^ 1 (S* (T A B)°). We will 
denote by P and Q the complementary projectors defined by this decomposition, that is, 

P b : % A A -» T h A B and Q b : T h A A -» F h , for all b e B. 

Theorem IV. 5. Let (L, B) be a regular constrained Lagrangian system and Rl be the solution 
of the free dynamics, i.e., iR L £lh = and in L 4>o = 1. Then, the solution of the constrained 
dynamics is the SODE R n h obtained by projection R n h = P{Rl\b)- 

Proof. Indeed, if we write R n h{b) = J?i(b) — <5b(-Ri(b)), for b G B, then we have 

*iWb)^L(b) = i RL{h) n L (b) - i Qb{RL{h)) n L (b) = -i Qb(RL{h)) n L (b) 

which is an element oiS*{T A B)°, because Q h (R L (b)) £ F h and <£o(Qb (#£(*>))) = 0. Moreover, 
using this last fact and since Rl is a SODE and S(Q(Rl)) = 0, we have that R n h is also a 
SODE. □ 

Let (x l ) be local coordinates on an open subset U of M and {eo, e a } be a local basis of sections 
of the vector bundle t^ 1 (U) — ► U adapted to 1^. Denote by (x l , y°, y a ) the corresponding local 

coordinates on A and suppose that the equations defining the constrained subbundle B as an 
affinc subbundle of A are 



* a = MS + = 0, for a e {!,..., r}. 
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Then, the local expression of the projector over T A B is 

P h = Id- J2 C ac (b)Z c (b)<E>(d TAA ^ a )(b), 

l<a,c<r 

for all b S B, (C a b) being the inverse matrix of (C ab ). 

Hence, if the constrained Lagrangian system is regular, the solution of the constrained dynamics 
is 

Rnh = R L \ B - C ac p T £(R L )(* a )Z clB . 

l<a.c<r 

From the regularity of the local matrices C we deduce that (P, Q) may be extended (in many 
ways) to an open neighborhood of B. Therefore, R n h may also be extended to an open neigh- 
borhood of B. This fact will be used in the following proposition. 

Proposition IV. 6. Let (L, B) be a regular constrained Lagrangian system, Ol be the Poincare- 
Cartan 1-section and C T A be the Lie derivate on T A A. 
i) If R n h is the solution of the constrained dynamics 



ii) We have that 



L R nh ®L-d L — L q{Rl) Vi 



Proof, (i) It follows since R nh = P{R L ) = R L - Q{R L ) and C^ L A L = d T A L. 

(it) Since Q{R L ) = ELi A ^a, with A Q = C ac p T £(R L )^ c ), and using Q3J) and we 

deduce that 

r 

c l*R L )®L=^^z a d TAA ®L+d TAA i^ KaZ Q L = -Y,^aiz a tt L = -A af i a J a . 



a=l 



Thus, (c rA Q L ) w e r(Tw^o). □ 



IV. 1.2. Cosymplectic Projection to H . We have seen that the regularity condition for the cons- 
trained system (L, B) can be equivalently expressed by requiring that the subbundle H is a 
cosymplectic subbundle of (T A A, SIl, <fio). Thus, the restriction (Q^,(Pq) to H of the cosym- 
plectic structure (£Il,4>o) is a cosymplectic structure on H. Therefore, there exists a unique 
solution on H of the equations 

i x n% = and i X (j>" = 1. (48) 

Proposition IV. 7. The solution of equations |^<S[ ) is precisely the solution of the constrained 
dynamics R n h- 

Proof. Since R n h is the solution of the constrained dynamics, then R n h(b) 6 T^B, for all 
b £ B. Moreover, using that R n h is a SODE, we obtain that R n h(b) 6 Gb- Then, R n h (b) £ i?b 
and it is obvious that it verifies 1®. Note that i Rnh Sl L (b) E S*{T^B)° = C H£, for all 
beB. □ 
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On the other hand, we have a direct sum decomposition 

T A A = H h ®H^, 

for all b € 6. We will denote by V and Q the complementary projectors defined by this 
decomposition, that is, 

V h : T h A A -> H h and Q b : T h A A -> H£, for all b e B. 

Then, we have the following result. 

Theorem IV. 8. Let (L, B) be a regular constrained Lagrangian system and Rl be the solution 
of the free dynamics, i.e., — and in L 4>o = 1. Then, the solution of the constrained 

dynamics is the SODE R n h obtained by projection R n h = V(Rl\b)- 

Proof. If b G B and X G T^A then, using that Qb(R L (b)) G H^, it follows that 

(ip b (ii £ (b))nf(b) + ^(b)(p b (iJL(b)))^ r (b))(p b W) 

= -(is b (ii,(b))^(b) + ^o(b)(Qb(i2 i (b)))^(b))(P b (X)) + ^(b)(P b (X)) (49) 

= ^(b)(V h (X)), 
Thus, we deduce that 

In particular, from Proposition IIV. 7l we obtain that 

1 = 4>"{Rnh) = ^(T(R LlB )). (50) 
Therefore, using (|49|) . we conclude that 

i^^nf = 0. (51) 

Consequently, from ([50]) and (fS"Tj) . it follows that 

□ 

Next, we will denote by the Poisson bracket on A induced by the algebraic Poisson 

structure Kj, given by 

for G C°°(A). 

On the other hand, we introduce the Hamiltonian section X^ a G T{t~£) associated with the 
function ^> a with respect to the cosymplectic structure (ClL,<f>o) which is defined by 

= \>- L \d TAA V a ) - (<F* A * a )(R L )R L , foraG{l,...,r}, 

where \>l ■ T A A — > (T A A)* is the vector bundle isomorphism defined in Note that, from 

(|40| and ijgjl . we obtain that 
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Then, using that {(grad ^ a )\ Bl Z a = —S(grad 1 i' a )\ B } a= i t ,,^ r is a local basis of sections of H- 1 , 
we deduce that the local expression of the projector V is 

V = Id- C ab C cd ([^\^ d } L + {d TAA ^ b ){R L )(d TAA ^ d )(R L )^ 

l<a,b,c,d<r 

+ J2 C ab (x^ a + (d T ^ a )(R L )R L )®S*(d T ^ b )- J2 C ab Z a ®(d TAA V b ), 

l<a,b<r l<a,6<r 

along the points of B. 

IV. 1.3. Poisson Projection to H. Assuming that the constrained Lagrangian system is regular, 
we have a direct sum decomposition 

T h A A = H h (BH^' AL , for all beB, 

where we recall that H- L ' Al = §a l (H°). We will denote by V and Q the complementary 
projectors defined by this decomposition, that is, 

V b : T A A -» H b and Q b : T A A -► H^ Al , for all b e B. 

Theorem IV. 9. Let (L,B) be a regular constrained Lagrangian system and Rl be the solution 
of the free dynamics. Then, the solution of the constrained dynamics is the SODE R n h obtained 
by projection R nh = V(R L \ B ). 

Proof. Suppose that b e B. Since Q b (R L (b)) = Q b (R L (b)) E F b = = G^' Al C h£ Al , 
V h (R L (b)) £ H h and R L (b) = V b (R L (b)) + Q b (R L (b)), we conclude that 

P b (R L (h)) = V b (R L (b)) = V b (R L (b)) = Rnh(b), 
Qb(R L (b)) = Q b (R L (b)) = Q b (R L (b)). 

□ 

Since also {Z a ,X^a} a =i,...,r is a local basis of H ] -' Al , we obtain that the local expression of 
the projector V is 

V = Id- J2 C ab C cd {* b ,* d } L Z a ®S*(d TAA * c ) 

l<a.b,c,d<r 

+ J2 c ab x A i®s*{d TXA ^> b )- J2 c ab z a ®{d TAA ^% 

l<a,6<r l<a,6<r 

along the points of B. 

IV. 2. The constrained Poincare-Cartan 2-section. Let (L,B) be a regular constrained 
Lagrangian system on a Lie affgebroid A of rank n with regular Lagrangian function L : A — > M 
and with constraint subbundle B of corank r. The equations for the Lagrange-d'Alembert sec- 
tion R n h = P(Rl\b) can be entirely written in terms of objects in the Lie algebroid prolongation 
(t^ 8 : T A B -» B, [; -1 T ~,Pj) of the Lie algebroid (tj : A -> M, [•, -]j,Pj) over the fibration 
Tg : B — > M. In order to do this, for every point b e B, we define 

w(b) = n L (b) - (i Qb(KL (b))^L(b)) A 0o (b). 
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uj is a section of the vector bundle A 2 (T A A)\ B — > S. We also have that <fo(b) A uj n (b) ^ 0, for 
all b e B. Thus, there exists a unique section X of T a A\q — > £> such that 

= and ix</>o = 1- (52) 

Note that X is the solution of the constrained dynamics. In fact, since (f>o(P(RL\B)) = 1) it 
follows that 

ip(R LlB )U = ip(B, LlB )^L + iQ(Ft, LlB )tt L - (iQ(R LlB )^L)(P(RL\B))<Po\B 
= -(*QCRi,| B )^l)0 P R £|b)) < A)|8. 

Thus, using that Q(Rl\b) is a section of F — > B and the fact that P(Rl\b) is a SODE, wc 
conclude that ip(R LlB )UJ = which proves that X — P(Rl\b) — Rnh- 
Next, we get the following. 

Theorem IV. 10. If Co and </>o are the restrictions of lo and 4>o to the vector subbundle T A B 
ofT^A, then the solution R n h = P(Rl\b) °f the constrained dynamics verifies the equations 

ix& — and ix<fto = 1- (53) 

Moreover, the unique SODE X on T A B satisfying (5$ is just R nh = P(R L \b)- 

Proof. Since the section R n h = P(Rl\q) satisfies ([52]) . then it also verifies (f53| . 

Now, let X be a SODE on T A B such that ixCj = and ix<j>o — !• Then, we have that 

(i x to)(P(Y)) = 0, for all Y section of T A A\ B -> B. (54) 

On the other hand, using that S(Q(Y)) = 0, <f> (Q(Y)) = and the fact that X is a SODE, wc 
obtain 

(i x oj)(Q(Y)) = -(i Q(Y) n L )(x) - (i Q{RL) n L )(x)MQ(Y)) + (i Q[RL) n L )(Q(Y)) = o, (55) 

for all Y section of T A A\ B -> B. Finally, from dM]) and IJ55J), we conclude that ix^> — which 
implies that X = R n h- □ 

Definition IV. 11. The 2-section Co is said to be the constrained Poincare-Cartan 2-section. 

Remark IV. 12. Note that (u),<f)o) is not a cosymplectic structure on T A B so that it may be 
another solution of the equations 

ix& — and ix<t>o = 1- 

V. Constrained Hamiltonian Systems and the nonholonomic bracket 

V.I. Constrained Hamiltonian Systems. We now pass to the Hamiltonian description of 
the nonholonomic system on a Lie affgebroid. The equivalence between the Lagrangian and 
the Hamiltonian description of an unconstrained system was discussed in Section [lI.61 We now 
consider a nonholonomic system, described by a hyperregular Lagrangian L and a constrained 
afHne subbundle t b : B —> M of the bundle T4 : A — * M. Denote by B the image of B 
under the Legendre transformation, which is a submanifold of V* . If B is again locally defined 
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by r independent functions \& a , then the constraint functions on V* describing B become 
iP a = ^ a o leg^ 1 , i.e., 

where the local expression of the Hamiltonian section is h(x l , y a ) — {x l , —H(x 3 , yp), y a ). 
A curve 7 : / — » V* is a solution of the equations of motion for the nonholonomic Hamiltonian 
system (h,B) if legj^ 7 : I — > .A is a solution of the Lagrange-d'Alambert equations for the 
system (L,B). Thus, using ((501) . (EU an d (|55|l. we deduce that 

7:/^*, t~7(*) = (*'(*),*«(*))> 

is a solution of the equations of motion if 

dx l i i dH 

-dT = Po + P "dy- a ' 

^ + i^L + ((T + —)v - —A H at3dr 
dt +Pa dxt +{Ca0+C <*Pdy^- XaH dyp' 

V a (^',y a ) = o, 

with i G {1, . . . , m}, a G {1, . . . , n} and a G {1, . . . , r} and where A^ are the Lagrange multipliers 

/ d^H \ 

and ?i a P are the components of the inverse matrix of the regular matrix (Tt a g) = ( — ^ — )• 

\dy a dy p J 

Note that 

(J?-H a/S )(xi,y 7 ) = (—oleg^ix^y,) = o leg^x^y,). (56) 

An intrinsic description of the equations of motion is obtained as follows. Consider the sub- 
bundle G C T^V* B -» B such that G 5 = {Tleg L ){G h ), with b = leg L (h) and b G B. It is easy 
to prove that is locally generated by the independent sections 

H a = ^H^^-^e°), for all a 6 {!,..., r}. (57) 
dy a dyp 

Note that, from (gTJ) and (O, it follows that 

(Tleg L ,leg L y({l a ) = fi a a 9 a , for all a. 

Moreover, it is clear that the restriction : B — > M to S of the vector bundle projection 
Ty : V* — > M is a fibration. Thus, we can consider the Lie algebroid t~? : T A B — » B and the 
Hamilton equations of motion of the nonholonomic system can then be rewritten in intrinsic 
form as 

= 1, (58) 

(rift., 77) being the cosymplectic structure on T-^F* defined in (flU]) and (jT5J) and rgo being the 
vector bundle projection of G° — > i3. It is said that the constrained Hamiltonian system (ft,, B) 
is regular if the Hamilton equations have a unique solution, which we will denote by R n h- 
Now, denote by b/, : T" 4 !^* — > (T* 4 "^*)* the vector bundle isomorphism given by 

K(x) = i x n h + v(x)v, for X G r x v*, 
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and by the (algebraic) Poisson structure on T A V* defined by 

A h (a,f3) = n h (b- 1 (a), ^ 1 (/?)), for a,[3 G (T"*V*)*. 

We also can transport the vector subbundle F C T A A\s — > £> and obtain a subbundle C 
T A V?g — > B such that F 5 = (Tleg L )(F h ), with b = Zegi(b) and b e 6. Moreover, it satisfies 
that 

4 = ^^)=^=^, 

where G 1 " (respectively, G ± ' Ah ) denotes the orthogonal to G with respect to the cosymplectic 
(respectively, Poisson) structure (flhtV) (respectively, Ah). Notice that F is locally generated 
by the sections Z a — X A a G r(r^ s ) characterized by the conditions 

%2 a ^lh = fl a and igT] = 0, 

for all a G {l,...,r}. Of course, Z Q and Z a are (TlegL, ^e<7i)-related. Moreover, its local 
expression is the following 

We will denote by (C ab ) the matrix which elements are 



Z a = -U a ^e a . (59) 



r b = (d T * v 'r)(z b ) = ~Ax^) = -H*??f-%£, on) 

oyp oy a 



for all a,b e {l,...,r}, where -X^J is the Hamiltonian section of ip a with respect to the 
(algebraic) Poisson structure A h , that is, X A £ = -jj A „ (d TAy * V a ), tU h : (T-*V*)* -> T'V* 
being the vector bundle morphism given by §A h (a) = i a ^h, for a G (T^V*)*. It is easy to 
check that C ab o Zeg L = C ab . 

In a similar way that in the Lagrangian side, we can consider the subbundle H C T A V*^ — > B 

such that -ffg = TP^SnGg = (TlegL){H\,), with b = ZegL(b) and b G £>, and its orthogonal J5" 1 - 
(respectively, ff- L ' Ah ) with respect to the cosymplectic (respectively, Poisson) structure (Qh,v) 
(respectively, A^). 

Using Theorem IIV.3I and the fact that the Lagrangian L is hyperregular, we deduce the next 
result. 

Theorem V.l. The following properties are equivalent: 

(1) The constrained Lagrangian system (L, B) is regular, 

(2) the constrained Hamiltonian system (h, B) is regular, 

(3) the matrices C = (C ab ) are non-singular, 

(4) T£B n F 5 = {0}, for all beB, 

(5) H£ n fi£ — {0}, for all beB, 

(6) H h n H^ Ah = {0}, for all b G 23. 

Assuming the regularity of the constrained system, we have the corresponding direct sum 
decompositions of T a V*\q (similar to the Lagrangian case). But we will only develop the 
Poisson decomposition. Since the condition (6) of the above theorem, we have the decomposition 

TgV* = ff s © H^ Ah , for all beB, 
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and, hence, one can define two complementary projectors 

: T^V* -» ff s and Q s : T^V* -> H^ Ah , for all beS. 

Moreover, (V, Q) is (Tlegz, ie<7i)-related with (P, Q). Then, the section V{Rh\g) is the unique 
solution of the constrained Hamilton equations (f55|) , where Rh is the solution of the Hamilton 
equations for the free dynamics, that is, R n h = F(Rh\s)- 

Now, we are going to construct a decomposition of T A A + , which allows us to obtain the solution 
of the constrained Hamiltonian system as the projection by T/x of a certain section. Given a 
Hamiltonian section h : V* — > A + , one can construct an affine function Fh ■ A + — > M with 
respect to the AV-bundle fi : A + — > U* as follows. For each ^ G „4+, with a; £ M, exists a 
unique Fh(<f x ) & R such that 

h(n((p x )) - fx = F h {(px)^A{x). (61) 

The function i*X : _4. + — * K is locally given by 

F h (x\y ,y a ) = -H(x\y a ) -y . (62) 

Note that this function was introduced in another way in Section III. 51 Using ^ , (fT"3"| , (fBl 
and (fUTj) . it is easy to prove that 

where is the Liouville section associated with the Lie algebroid A and (7~/i, /j.) is the Lie 
algebroid epimorphism between the Lie algebroids r ^ A+ : T A A + — > .4 + and r~ / : T-^U* — > U* 

defined by T/ifS,^.) = (5, (T^ps^ J), for (5,X V J e T A A+. As consequence, we have 
that 

n J =(T f M,n)*Q h + 9 h A(T[M,n)*ri, (63) 

where 6h = d rAA+ Fy l and f2 j is the canonical symplectic section associated with A. 
We will denote by E S T(t^ 4+ ) the section given by 

E(<p x ) = (0(x),-l%(<p x ))eT A A + , for all ^6.4+ (64) 
where 1^ G X(.4 + ) is the vertical lift of the section I.4 e T(t4+). Using (|62|) . we obtain that 

kerT Vx ti=<E(<p x ) >, 6 h (y x ){E(y x )) = 1, (65) 
for all (/3a; € A^T, with a; S Af. 

Now, we introduce the vector subbundle (T A V*) H C T w4 ^l + — » .4+ whose fibre at point 
(p x e At is 

(^i.)^*)^. = e = 0}. 

It is not difficult to prove that [T Vlc y) ur A V , )H ■ (X^L )U*)^ — * T-4 ,y* is a linear 
isomorphism, for all tp x £ A^ ■ Then, we deduce that 

T*A+ = (Tfi Vm) vX.® < E(<p x ) >, (66) 
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for all ip x e A+, x G M. We will denote by ^ : , V* -> {T$ Vx )V*)% x the inverse map of 

v«*ti\{T$^v)* m ■ vi*.) v X. -> T ^) v *- So ' if x e !t follows that 

X = (T V ^)(X V X, + 9 h {<p x ){X Vm )E{v 9 ). (67) 
On the other hand, consider the section fj G r((r~* + )*) defined by 

f\{<f> x ){a,X,p x ) = 1.4(a), for ^ G .4+ and (a,X p J G 
A direct computation, using @ and |(55]). proves that 

a(x%) = -i, 

where -X^jf G T(t~ 4 ) is the Hamiltonian section of Fh with respect to the symplectic structure 
f2jj. Thus, since (T(J,,{i)*r) — fj, we deduce that 

)? (M^))((r^ Ai )(^(^))) = -i. 

Therefore, from (|6"5|) . we obtain that 

(i (7;; ^ )(x ^ ( ^ )) 0U//(^)))((r^ M )(^J) - 0, for Z Vm G T^+. 

This implies that 

{%^){X^{ip x )) = -RhMtp*)), for <p x eA+. (68) 

Now, if a G r((T^ v ')*), then the Hamiltonian section X^ h G ^( T j^) °f a with respect to the 
(algebraic) Poisson structure is characterized by the following conditions 

i x A h Q h — a - a(Rh)r) and i x A h T) = 0. 

On the other hand, denote by G T(A 2 (r^ 4+ )) the (algebraic) Poisson structure associated 
with the symplectic structure flj. In other words, 

A x (a, 0) = (a), b- 1 08)), for a, G r((rj + )*), 

b : T-^^l + — > (T-^^l + )* being the vector bundle isomorphism induced by f^. Note that if 
<5 G r((Tv> + )*) then 

= -Ka^(5). 

Proposition V.2. TTie prolongation of fi 

2> : T X A + -> TV* 
is a Poisson morphism over fi, that is, if a G r((r~ A )*) i/ien 

T M°^ ( n / M , M) . a = ^°M. (69) 
Proof. If (jOj; G ,4+ then a direct computation, using (j4]) and l|14jl . proves that 

»?(M(^))[(^.A*)(^ M) . (*'x))] = o. 
Thus, from (J63j> and ((68]), it follows that 

(i , ^fi/i(M(^)))((^>xM)(^*)) = a OOz))((^*M)(^>J) 
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+0^)(^,^ ) ,j¥'x),^(^))r?(M(^))((^M)(^J) 

for Z Vx G %£A + . This implies that 

V u , „^/.(m(Vx)) =a(M(^x)) -a(n((p x ))(R h (fx((p x )))r](fi((p x )). 

Therefore, we deduce that ([MI) holds. □ 
Remark V.3. From Proposition IV. 21 it follows that 

Next, we consider the submanifold B + of A + defined by B + = and the vector subbundle 

H+ C T J A~l B + -> 23+ given by 

^ = (^-m)~ 1 (^.)) = (^(v-))^®< £ (^)> 1 fora11 V»eB+ (71) 

where = CW)" 1 n (^V*)^. 

We will denote by (i/+J ±,n -4 the orthogonal to if+ with respect to the symplectic structure 
n A (ip x ) of T*A + , for all tp x G B+. In other words, 

(H$.) ±>a * = l-HW) = U A ((H+J). (72) 

Using ([7T|) . we have that 

)° = (T^ M r((;w))°) 

and, from ((TDJ and ([72]). it follows that 

(T^/x)((if+ = (i?^)^, for all p x G B+ (73) 

In addition, using that {/2°,c2 r is a local basis of H°, we obtain that {(T/z, fJ-)*n a , 

(Tu,/i)*(d r ^ y > Q )}isalocalbasisof (H+)°. This implies that {X^ w - a ,X^ , } 

is a local basis of (H + ) ± ' n ^. 

Note that, from Proposition IV. 21 we have that 

(T^iX^^-M) = Xp{^ x )) = Z a (^ x )), 

< T ^><4^^W^ M = X £WP.)). (?4) 
Then, we will prove the following result. 

Theorem V.4. The following conditions are equivalent: 

(1) The constrained Hamiltonian system (h, B) is regular, 

(2) i?+ n (H+ m ) ± > a x = {0}, for all ip x G B+ , with x G M. 



Proof. (1) (2) Suppose that the constrained Hamiltonian system is regular. Then, since 
B + = /i _1 (z?) and using Theorem IV. 1[ we have that n H^ v h ^ = {0}, for all cp x G 

B+. If X G H+ m n (if+ then, from ([711), * = (T Vm /i){X) G i^( Vx) . Moreover, using 

([751) . we deduce that X G ( J ff A1 ( Vx) )- L ' A ' 1 and, therefore, X = 0. Thus, X G < E(y> a ) >. 
Therefore, from ([Ml and ([MI), ixftjffax) = A((T //, fx)*r])(cp x ) G (#<+ )°, with A G E. However, 



34 



D. IGLESIAS, J. C. MARRERO, D. MARTfN DE DIEGO, AND D. SOSA 



((Tfj.,fM)*rj)(<f x ) g {H+J°. In fact, there exists R n h(fx) £ H+ x such that (Rnh(<Px)) = 

Rnh(p(<Px)) and, thus, ((Tfi,(i)*rj)(tp x )(Rnh((Px)) = il(tJ'( ( Px))(Rnh(fJ>(<Px))) = !• Consequently, 
A = and X = 0. 

(2) => (1) Suppose that H+ x n (H+J^X = {0}, for all ip x £ B+ and take X £ nH^S. 
Using [[75]) . there exists X G such that {T Vic ^){X) = X and, since X £ H^ v ^, we 

have that X £ H^. Thus, X — 0, which implies that X = 0. Finally, using Theorem IV. 1( we 
conclude the result. 

□ 

Assuming that the constrained system is regular, the bundle H + — > B + is a symplectic sub- 
bundle of the symplectic bundle (T a A + \b+ — > £> + , (fi j)|b+) and we have a direct sum decom- 
position 

T*A + = 8 (il+J-^, for all ^6B+ with i £ M. 
Let us denote by V + and Q + the complementary projectors defined by this decomposition, that 

7>+ : 7^.4+ -» ff+ and Q+ : 7^.4+ -» ) x > n A, for all <p x eB+. 

Theorem V.5. Let (h, B) be a regular constrained Hamiltonian system. Then, the solution of 
the constrained dynamics is the section R n u obtained as follows R n h o fj, = — 7~(j,(V + (Xp A )). 

Proof. If f x G £>+ then, from (TTTj) and since the map 7^ x /i : T^A + — > 7j5 ^V* is a linear 
epimorphism, it follows that 

{T Vm ii)(H+ 9 ) = H Kvm) . (75) 
Thus, if Z Va; £ Z£A+ then, using (|73[) and (|75p . we deduce that 

C?;,m)(^(^J) e h^.), (r^ M )(QX(^)) e ( j5 M*>.)) x ' Ah 

and it is clear that 

CWX^J - (^)(7>+ (Z v .)) + (^/i)(Q+ (V)). 
Therefore, we have proved that 

i.e., 

° Vt» = ^(^) ° T ^ and r f^ ° 2jx = CW*) ° V.M- ( 76 ) 
Finally, from |68|) and (|76|) . we conclude that 

Rnh(KVx)) = ^)(i2fc(M(Vx))) = -^(V^W)) = -^.^.(^(Vx)))- 

□ 

Next, we will denote by the Poisson bracket on V* induced by the algebraic Poisson 

structure given by 

tyrfh = -{d TXv *i>'){X^) = Sl h (X$»,Xfr), (77) 
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for ip,ip' G C°°(V*). 

Then, using ((HO)). (|7T|l . dJ) and the fact that /2 a (Z 6 ) = (see (H7J and JSS])), we deduce that 
the local expression of the projector V + is 

v+ = / c ^-(c ah o M )(c cd o M )({^^^} ^l o M )x^ ) ,_ o 0(r M , Ai )^ c 

(C ab o M )X ( % M) ,_ a ® (r Ml/ i)*(cF^> b ) (78) 
+ (C a6 o M )X n ^ ®(T M , M )*/i 6 , 
along the points of B + . 

Now, we will denote by Oh the section of (T A V*)*^ — > B given by 

^(b) = -fe nh ^)(b), for all be B. 

Then, using ((711), Theorem EH and the facts that Im(V + ) C Im{Q+) C (iJ+) ± ' f2 ^, 
■Rnft G -ff and ijy h T) — 1 , it is not difficult to prove that 

^(m(^))(Sm^)(^m(^J)) - -Mv*)(Q£. (79) 

for all X Vx G T A x A + and ip x G B+, with a; 6 M. Note that r){ii{ip x )) o 7^.^ o = n(fi(tp x )) o 

Qn(ip x ) ° — 0- We also recall that 6^ = d TAA+ F h . 
Thus, using l[6T|). ([76]) and ((75)1. we obtain that 

for all X Va; e T^„4+, with ^ S B+ and a; G M. 

V.2. The nonholonomic bracket. We consider a regular nonholonomic system on a Lie aff- 
gebroid A described by a hyperregular Lagrangian function L : A — > R and a constraint affine 
subbundle Tjg : B —> M of the bundle T4 : A — > M. We will denote by (7i, B) the corresponding 
regular constrained Hamiltonian system, by h the restriction to B of h, by B + the submanifold 
of A + given by B + = ^ _1 (B) and by 

: - 21+ and Q+ : T*,4+ - (22+ )^~, 

the corresponding complementary projectors, for 99^ G B+. 

We have that B + is the total space of an AV-bundle over B. In fact, the affine bundle projection 

is the restriction : B + — > B to B + of the canonical projection // : A + — » V*. 

Note that the restriction rg+ : B + — ► M of the projection T4+ : ,A + — > M is a fibration. Thus, 

one may consider the prolongation T A B + of the Lie algebroid A over t b + : B + — * M. T A B + 

is a Lie algebroid over B + . 

Now, suppose that 

h\h" : B ^ B+ eT{fi B+ ) 
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are two sections of the AV-bundlc yu e + : B + — > B and that 

h',h" :V* ^ A + eT(ij>) 

are arbitrary extensions to V* of h' and h", respectively. 
We will denote by 

F h ,,F h „ : A + — > R 

the affine functions associated with the sections /i' and h", respectively, and by X F £ and X F * t 
the corresponding Hamiltonian sections. Then, we have that 

h'(fj,(ip x )) - (fix = F h >(tp x )l A (x), h"(fj,(ip x )) - tp x = F h ,>(tp x )l A (x), 

Q A = F h' > i x n A Q A = dTAA+ F h" , 

F h> F h" 

for f x e At- 

Moreover, we will prove the following results. 
Lemma V.6. The real function on B + 

does not depend on the chosen extensions h' and h" ofh! and h" , respectively. 
Proof. If <p x e B+ then Q+ (X^(<p x )) e (H+J^a and, thus, 

^x)(Q+ {x n F *{y x )),v+ x = 

which implies that 

^ + (<1, B+ )^ + (^| B+ )) = ( 81 ) 
Now, let ft^ : V* —* .4+ be another extension of ft,' : B — > £> + . Then, it is clear that 

- ^/»' 1 )|B+ = 0- 

Therefore, since V+ x (X F *(ip x )) € C T^B+, we obtain that 

^ x )(x^> x ) - fa*))) = p£ + (^)))(^' - fy) = °- 

This proves that 

«# + (^ |B+ )^ + (^ |B+ )) = ^ + (<t |B+ )^ + (^ |B+ ))- 
Consequently, using that is skew-symmetric, we deduce the result. □ 
Lemma V.7. The real function on B + 

is basic with respect to the affine bundle projection : B + — ► B. 
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Proof. Denote by 9h> the section of (T A A + )* — ► A + denned by 

9 h , = <f^ + F h , 

Then, from ([81]), it follows that 

^(V + (X% lB + ),V+(X%, {B +) ) = 6 h , ]B+ (X^ r , {B + ) - O h/lB+ (Q + (X^ iiB+) ). 

Moreover, if Il_4+ is the linear Poisson structure on A + induced by the Lie algebroid A, we 
have that (see ^) 9^ (X F * rl ) = {Fh> , F^n }n A+ ■ In addition, the vertical bundle of the fibration 
|U : A + — > V* is generated by the vertical lift 1^ of the section 1^ G r(r4+). Since is a 
1-cocycle of A, we deduce that 1^ is an infinitesimal Poisson automorphism of II_4+ and 

l%({Fh>,Fh»}n A+ ) = {-^(F h '),F h ,,}n A+ +{F h ,AV(F h ,,)} nA+ 

= -{l,F h „}u A+ -{F h ,,l}u A+ =0. 

This implies that the function 9h>(X F ^ ) is basic with respect to the fibration /i : A + — * V™. 

Thus, the function 9h'{X F ^ ii )\Q+ is basic with respect to the fibration : B + — > 23. 
On the other hand, from (|80p . we deduce that 

M*>x)(Q+ (*{£>»))) = ^'(¥'x)((^)(T ¥; ^(X^(^))))^) 
for <p x G B+. 

Now, using (foTj) and one proves that 

Mv«)(£fe>*)) = = -4 0M¥>x) -*) = !, for ^ G At- (82) 

at |t=o 

Furthermore, proceeding as in Section rV.il (see (|68j)). we have that there exists a section G 
Tirf) such that 

{%p m it){X%* { Vx }) = -R h „{»(<p x )), for all <p x G A+. 

Therefore, 

M^(<& = -^(^)((4(^)(^"(m(^))))?J 

+Oh{^x)){Q^){Rh"{^x)))), for ^ G S+. 

Consequently, it only remains to prove that the function 9f l '(Q(Rh")) H ■ B + — * K- defined by 

9 h ,{&{R h »)) H {v x ) = e h ,{y x ){{^ {Vx) {R h »{»{v*))))Z x ) 
is /xg+-projectable or, equivalently, that 

C T E *f(9 h ,(C>(R h „)) H )=Q. 

In fact, we will see that if Z G F(t^) and Z ff is the section of T A A + given by Z H (tp x ) = 
(Z(M(¥>x)))£, then 

£^- 4+ (^(Z H )) = 0. (83) 



38 



D. IGLESIAS, J. C. MARRERO, D. MARTfN DE DIEGO, AND D. SOSA 



Note that, from {82]), it follows that 

£ T / A+ 9 h , = <f A * + {cf* + F h ,) = d^ + {6 h ,{E)) = 0. 

Thus, 

ct A+ {e h ,{z H )) = e h ,mz H Yf). 

Now, since d rAA+ 9h — and using that 9h(E) = 1 and the fact that 9h(Z H ) = 0, we obtain 
that 

8 h {\E,Z H ff)=0. (84) 

In addition, 

T/i o E = 0, T^l o Z H = Z o fj, 

and the pair (T/i, u) is a morphism between the Lie algebroids T A A + — > A + and T A V* —* V* . 
This implies that 

Tao{E,Z H ff =0. (85) 
Therefore, from {67]), (J84l> and {85]), we deduce that [E, Z H f~£ + = and {83]) holds. 

□ 

From Lemmas IV. 61 and I V. 7( we have that there exists a real function {ft/, h"} n h £ C°°(B) which 
is characterized by the following condition 

{l',ft"Uo fr ,!i I (P+(I^ |8+ ),P+(X^ |B+ )). (86) 

This function is called the nonholonomic bracket of the sections ft' and h" and the 

resultant map 

{vU:rWxrM^c»(8) 

is called the nonholonomic bracket associated with the regular constrained Hamil- 
tonian system (h,B). 

Theorem V.8. a) The nonholonomic bracket {•, •}„/,. : r(/i e +) x r(/ig+) — > C°°{JS) associated 
with the system (ft, 6) is an almost-aff-Poisson bracket on the AV-bundle u&+ : £> + — > B, that 
is, 

(i) {•) '}nh is a bi-affine map, 

(ii) {•, - } n h is skew- symmetric and 

(iii) Ifh! G F( ( u e+ ) iften 

{ft', -}„ ft : T(M5+) -» C°°(i3), ft" ~ {ft', ft"} nh , 
is an affine derivation. 

b) If f G C°°(/B) is an observable and {h, : C°°(23) — ■> C°°(B) is the linear map associated 
with the affine map {ft, -} n /j : F(/ig+) — » G°°(p), then 

f={hj} a n l h 

where f is the evolution of f along the solutions of the constrained Hamilton equations. 
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Proof, a) It is clear that {•, - } n h is skew-symmetric. 

Next, we will prove that {■, -}„h is a bi-affme map. Suppose that ft' and ft" are sections of the 
AV-bundle : B + — > B and that /' and /" are real functions on B. If ft', ft" E are 
arbitrary extensions of ft' and ft,", respectively, and f',f" G C°°(V*) are arbitrary extensions 
of /' and /", respectively, then we can consider the real functions on B + given by 

^(<f |B+ ^ + (</c,) |B+ )) = -^( n /c,) |B+ ^ + (<1 |B+ )) 5 

^(/^) IB+ ^ + (^ IB+ )) = -^(<v^ +( </w )} ' 

Now, proceeding as in the proof of Lemma IV. 6[ we deduce that these functions don't depend 
on the chosen extensions ft', ft" and /', /" of ft', ft" and /', /", respectively. 
Thus, we can introduce the linear maps 



{h>, -ri ■ c°°{B) - c°°{b + ), /" - {h>J>%' h = ^(x^ |b+ ,p+(x^ om)|b+ )), 

{•^nk : c°°(^) - c°°{b+) 1 r - = nA) |Bt /(^, |Jt )) 

= -{ft"j%, 

and the bilinear map 

{vS, = C°°(£) x C°°(B) <7°°(B+) 

Using that 

%' + /«)=f/»" + (/"°M), 

we obtain that 

{ft', ft" + f"} nh o Mb+ = {ft', h"} nh o MB+ + ^Tf} a nh- 
This implies that there exists {h',f"}£ h G C°°(S) such that 

and, in addition, 

{ft', ft" + f"} nh = {ft', h"} nh + {ft', f"} a nl (87) 

Since {•, h"} l ° h = -{ft", -}° l h , we have that there exists {/', h"}% h G C°°{B) such that 
and 

{ft' + /', ft"} n/s = {h', h"}nh + {/', h"}nh- (88) 
On the other hand, from (JHZJ) and (jS5]>. it follows that 

{ft' + /', ft" + /""}„, o Mb+ = ({ft', ft"} nft + {ft', + {/', ft"}^) o ^ig+ + { 
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Therefore, there exists {f'J"}" h £ C°°(B) such that 

{f\n l Lo^B + ={f7f} l L 

and 

{h 1 + f, h" + ru = {h\ h"u + fh>, rri + {/', h"Y n a h + {/', 

This proves that the map {•, - } n h is bi-affine. 

Moreover, the linear map associated with the affine map {h',-} n h is just the map {ft',-}^ : 
C°°(B) -> C°°(B) and we have that 

for / £ C CC (B), where / £ C°°(U*) is an arbitrary extension of /. Thus, since the pair (Tfi, fi) 
is a Lie algebroid morphism, it follows that 

{h'jy&Mv*)) = -((r M , M )*(d r ^7))(^)^ + (^)(^)). (89) 

In addition, proceeding as in Section ED (see (|68|) ). we deduce that there exists a section 
R h > £ r(r~^) such that 

[T Va MX%£{<p x )) = -Rh.(jiM<P*)), fOT <P* G S+, (90) 
and, using (|76[) . we obtain that 

= -^ s+ (^)(^(MB+fe))) 5 for ^ G B+ (91) 
Therefore, from (f89l) . (f9Tjl and since T-^B is a Lie subalgebroid of T A V* , we conclude that 

=P^(fl/.'))(/). (92) 

Consequently, {ft', -}^ h is a vector field on B and {ft', -} n h is an afhne derivation. 

b) It follows using and the fact that the solution of the constrained Hamilton equations is 

the section V{Rh\g)- 

□ 

Now, let ft', ft" be two sections of the AV-bundle /J,g+ : B + — » B and ft', ft" £ T(fi) be two 
arbitrary extensions to V* of ft' and ft", respectively. Then, we may consider the nonholonomic 
bracket {ft', h"} n h of the sections ft' and ft" and the aff-Poisson bracket {ft', ft"} of the sections 
ft' and ft" defined in JT5J). Using ©, JI5J), id, ©, ([Ml) and (HOI), we obtain that {h',h"} nh 
and {ft', ft"} are locally related by the following condition 

{ft', h"} nh = ({ft', ft"} - C ah C c M\ i> d } h fl c (R hl> )fl a {R h> ) - C ah {4>\ h"} v fl a {R h> ) 

+C ab {iP a ,h'} v p, b (R h „))_, 

where {•, ft"}y '■ C°°(V*) — » C°°(V*) is the linear map associated with the affine map {•, ft"} : 
r(/i) — > C°°(U*) (see (j^U])). Thus, if the local expressions of the sections ft' and ft" are 
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h'(x\y a ) = (x\-H'(x : > ,yp),y a ) and h"(x\y a ) = {x\ -H"(x 3 , yp), y a ), then, using 
(fig]) . HU), Jill), |J22J), dST]) and (UU), we deduce that 

rT.rT.rn ^{H'-H") t .dH' dH" dH' dH" . ^ 7 d(H'-H") ^ dH' dH 

{h>, h»} nh = pi + ) + C 7 ^7 - - ^ 



+C a b 
—Cab 



dx* ^ ay dx l dy a dy a dx l ' Uq * 7 0j/ o a/37 %* <% 

^^_d^dj/_\ c7 ^ b djj c ^ d(H" - H) d^ a gv d{H' - H) 



Pa \dx i dy a dy a dx i J ' ~ a P a ~> dyp dy a J dy 1 '" dyp dy e '" dy u 
l di^_ j ( dip b dH" d^ b dH" \ d^ b ^ dH\-\dj) a ^ g d(H' - H) 

Po dx* Pa \dy a dxi dx * dya ) + dy a yj{ 0a a " dyp >\ 0y 7 dy p 



dip a i fd^ a dH' dip a dH'\ dip a dH' 

HaZa?" ~ ~dx T d^J + ~dZ y ^ 0a ~ a P~d~^' 



■H lfi - 

dx* t ' a \dy a dx* dx* dy a J ' dy a aiK ~ Ua ~ ati d Vfi '\ d Vl dy p 



di>* y0 d(H"-H) 



VI. Examples 

VI. 1. Lagrangian systems with linear nonholonomic constraints on a Lie algebroid. 

In this section, we will discuss the particular case of Lagrangian systems with linear nonholo- 
nomic constraints on a Lie algebroid (see [?])• 

First of all, we will recall an standard construction which will be useful in the sequel. 
Let te '■ E — > V be a Lie algebroid over a manifold V with Lie algebroid structure ([•, -Je, Pe)- 
Then, the vector bundle te ■ E x R — » V admits a natural Lie algebroid structure ([•, -]b x R) 
Pexr) given by 

[(V,/),(V, ff )] BxR = (lX,Y] E , PE (X)(g)-p E (Y)(f)), 
Pexr{X,J) = p E (X), 
for (V,/),(V,« ? )eF(r £ )xC 00 (V). 

Next, suppose that ry : V — > M is a Lie algebroid over a manifold M with Lie algebroid 
structure Then, the affine bundle T4 = ry : A = V —> M is a Lie affgebroid 

over M. In fact, the bidual bundle rr : A — ► M may be identified with the Lie algebroid 
fy : V x R — > M. In addition, it is easy to prove that the prolongation T A A of A over the 
fibration T4 : A — » M is isomorphic to the Lie algebroid fy v : T V V xl^y. 
We will denote by S : T v V -> T v V the vertical endomorphism on T V V (see [13 US]) and by 
E Q the section (0, 1) of fy v : T V V x R -> V. 

Now, let L : V -> R be a Lagrangian function on V and £ L = A(L) - L G C°°(V) be 
the Lagrangian energy (here, A is the Liouville vector field of V). Then, the dual bundle to 
t t £ : r^A -> A is isomorphic to the vector bundle {f v v )* : (T V V)* x R -> V and, under this 
identification, the 1-cocycle </>o is the section (0, 1) and 

il L =u; L + (0,1) A d TVy E L , (93) 

where u>l '■ V — > A 2 (T V V)* is the Poincare-Cartan 2-section associated with the Lagrangian L 
on V ([HH3]). 

Next, suppose that : U — > A'/ is a vector subbundle of ry : V — > M. Then, rrs — tjj : 
B = U — y M is an affine subbundle of T4 : .4 — > M and we can consider the affine constrained 
Lagrangian system (L, B) on ^4. Moreover, if X is a section of t- 4 = fy v : T A A = V xl^ 
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V then, from (|93p. it follows that X is a solution of the Lagrange-d'Alembert equations for the 
affine nonholonomic system (L, B) if and only if 

Y = Y\u = {X - E )\u is a SODE on U, 
{lyoJL - d rVv E L )\ u G T^.^v^o), (94) 

(|94[) are just the Lagrange-d'Alembert equations for the linear nonholonomic system (L, U) in 
the terminology of [7]. 

Thus, the affine constrained Lagrangian system (L, B) on A is regular if and only if the linear 
nonholonomic system (L, U) on V is regular in the sense of [7]. 

Therefore, using the results in Section ITVl of this paper, one directly deduces some results, for 
linear nonholonomic systems on Lie algebroids, which were obtained in [7j (see Section 3 in [7]). 
Now, assume that the Lagrangian function L : A = V — > R is hyperregular (that is, the 
Legendre transformation legL '■ A — V — > V* is a global diffeomorphism) and that the affine 
constrained Lagrangian system (L, B) is regular. Then, the space A + may be identified with 
the product manifold V* x R and, under this identification, we have that: 

(i) The vector bundle projection T4+ : A + — > M is the map f v : V* x R — > M given by 

fy (a K , i) = ry(a x ) = x, for (a x ,t) £ l^* x t and x G M. 

(ii) The fibration : A + —* V* is the canonical projection pr\ : V* x R — > V* on the first 
factor and, thus, the spaces r(/z) and C°°(V"*) are isomorphic. 

(iii) The Hamiltonian section h = LegL ° ^g^ 1 G T(fi) is the Hamiltonian energy H G 
C°°{V*) of the free system, that is, H = S L ° ^g^ 1 : V* -> R. 

(iv) The canonical aff-Poisson bracket on the trivial AV-bundle pr\ : V* x R — > is just 
the linear Poisson bracket on T^* induced by the Lie algebroid structure on V. 

Next, let B = U be the Hamiltonian constrained submanifold of V*, that is, B = U = legL(U). 
Then, the constrained AV-bundle : B + — * B may be identified with the trivial AV-bundle 
pri : U x R — * U. Therefore, the nonholonomic bracket associated with the nonholonomic 
system (L, B) is an almost-Poisson bracket on U, i.e., a R-bilinear map 

{;-}nh:C°°(U) x C°°(C7) ^ C°°(C7) 

which is skew-symmetric and a derivation in each argument with respect to the standard product 
of functions ({•, -} n h doesn't satisfy, in general, the Jacobi identity). 
Finally, since the restriction of the Legendre transformation to U 

legu = legh\u - U 

is a global diffeomorphism, one may consider the corresponding nonholonomic bracket on U, 
which we also denote by {•, - } n h, defined by 

{f,g}nh = {f °legu 1 ,goleg u 1 } nh olegu, for f,g E C°°(U). 

This bracket was introduced in [7] and its properties were discussed in this paper (see Section 
3.5 in 0). 
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VI. 2. Standard affine nonholonomic Lagrangian systems. Let r : M — > R be a fibration 
and A — J x t be the 1-jet bundle of local sections of r : M — > R. Then, as we know (see 
Section ril.3j) . the affine bundle ta = ti.o : A = J 1 t — > M is a Lie affgebroid modelled on the 
Lie algebroid ry = (7TM)|y r : V = ^ T — * M. Moreover, the bidual bundle A may be identified 
with the standard Lie algebroid ttm : TM — > M and, under this identification, the I-cocycle 
\a of A is just the closed 1-form r*(dt). 

Note that if (t, g a ) are local coordinates on M which are adapted to the fibration r then 
d d 

{eo = — , e Q = Tj - ^;} is a local basis of sections of the Lie algebroid = 7r^/ : .4 = TM — -> M 
such that lyt(eo) = 1 and l^(e Q ) = 0. Furthermore, we have that 

fj( e o) = -Q t ' P^( e a) = [eo, e a ];j = [e a , = 0, (95) 

for all a and /?. 

On the other hand, it is easy to prove that the prolongation r~* : T A A — > .4 of 4 over the 
fibration ta '■ A — > M is isomorphic to the standard Lie algebroid 7Tji T : T( J 1 r) — » JV. Thus, 
the space IXt- 4 ) may be identified with the set of vector fields on J 1 r. 

Now, let L : J 1 t — > R be a regular Lagrangian function. If ti = t o ti.o : J x t — » R is the 
canonical projection then, under the above identifications, the 1-cocycle <po of the Lie algebroid 
t~* : T-^4 — )► _A is the closed I-form rj — Tl(dt) on J 1 r. Moreover, the Poincare-Cartan 
2-section fi^ associated with L is just the Poincare-Cartan 2-form on J 1 t. 
Next, suppose that rg : B — > M is an affine subbundlc of T4 = n.o : 4 = JV — > M. Then, we 
can consider the affine constrained Lagrangian system (L,B). 

In addition, we have that a vector field X on J 1 r is a solution of the Lagrange-d'Alembert 
equations if and only if 

fof!t)| S er(r s , (Ti?) „), (W?)| B = 1, X| S e*(B). (96) 

These equations were considered in [T3] (see also [H]). Note that the first two equations imply 
that X is a SODE on B. 

From (l96|) . it follows that the affine nonholonomic Lagrangian system (L, B) on A is regular if 
and only if it is regular in the sense of [13]. Thus, using the results in Section [TVl one directly 
deduces some results, for standard affine nonholonomic Lagrangian systems which were obtained 
in [T3] (see also [TB]). 

On the other hand, we have that the space A + may be identified with the cotangent bundle 
T*M to M and, under this identification, the fibration /i : A + — > V* is just the dual map 
iy T : T*M — > F*r of the canonical inclusion iy T : W — > TM. 

Next, assume that the Lagrangian function L : A = J 1 t — > R is hyperregular and that the 
affine constrained Lagrangian system (L,B) is regular. Then, we can consider the Hamiltonian 
constrained submanifold B = legL(B) of V*t and the constrained AV-bundle fi B + : B + = 
(z* VT )-HB)^B. 

Now, denote by ft : V*t — > T*M the Hamiltonian section induced by the hyperregular La- 
grangian L (ft = LegL ° legj^ 1 ) and suppose that ft,' and ft" are two sections of the AV-bundle 
: 23+ — > B and that ft', ft" : V"*t — > T*M € r(«y r ) are arbitrary extensions of ft' and ft", 
respectively. Then, we will obtain the local expression of the nonholonomic bracket of ft' and 
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h". For this purpose, we take canonical local coordinates (t,q a ,p a ) (respectively, (t, q a ,p,p a )) 
on V*t (respectively, T*M) such that 

h(t,q a ,p a ) = {t,q a ,-H{t,q ,p ),p a ), 

h'{t,q a ,p a ) = {t,q a ,-H'(t,qP,p p ),p a ), 

h"(t,q a ,p a ) = (t,q a ,-H"(t,qf 3 ,pp),p a ) 

and the local equations defining B as a submanifold of V*t are 

1 p a (t,q a ,p a )=0, ae {f,...,r}. 

Then, using the general local expression of the nonholonomic bracket (see Section IV. 2j) and 



95|) . we deduce that 
{h',h"} nh 



rdH' dH" dH' dH" dH' dH" 
i~dt df + ~dq" dp a ~ ~dp~ a dq a 

C ah cJ—— - d ^ d ^ C -W d{ - H " ~ H) d ^ a H eu 9{H' - H) 
a c V da a dv n dv n da a J 9z>~ dvet dpa dv,. 



. dq a dp a dp a dq a J dp 1 dpp dp e dp v 

c dip b /dip b dH" dip b dH"\idip a n70 d(H' - H 
dt \dq a dp a dp a dq a )\ dp 1 dpp 



n -dip a (dip a dH' dip a 8H'^ «» Ab 

Cab 



\da a dvrv dvrv da a J J dn~ 



d^^ d(H"-H) 



dt \dq a dp a dp a dq a J J dp 1 dpp ) \b 

Finally, we remark that the linear-linear part of the bi-affme map {•, -} n h '■ r(/ig+) x r(/Xg+) — > 
C°° (8) is an almost-Poisson bracket on B. This bracket was considered in [3] . 

VI. 3. A homogeneous rolling ball without sliding on a rotating table with time- 
dependent angular velocity. [H [501 HH] Consider the vector bundle tjj : A — > M where 
.4 = TR 3 xl 3 , M = R 3 and tj is the composition of the projection over the first factor with 
the canonical projection 7r K 3 : TR 3 — > R 3 . Denote by (t,x,y) the coordinates on M. A global 
basis of sections of r j can be construct as follows 

( d d d ) 

where ui,u 2 ,u 3 : R 3 — > R 3 are the constant maps u\(t,x,y) = (1,0,0), u-2.it, x,y) — (0,1,0) 
and u 3 (t,x,y) = (0,0,1). Consider the projection pj : TR 3 xR 3 ^ TR 3 over the first factor 
and the Lie bracket on the space of sections T(t j) where the only non-zero Lie brackets are 

[e4,e 3 ]j = e 5 , le 5 ,e 4 jj = e 3 and |e 3 , e 5 jj = e 4 . 

Thus, ([•, •] j, p t) induces a Lie algebroid structure on A = TK 3 x R 3 . Denote by (t, x, y; i, x, y, 
u> x , LOy,oj z ) the coordinates on A induced by {eo, e\, e?, e 3 , e^, e$}. 

Moreover, 4> '■ TE. 3 x R 3 — > K given by <t>(t,x,y;i,x,y,uj x ,ujy,u; z ) = i is a 1-cocycle in the 
corresponding Lie algebroid cohomology and, then, it induces a Lie affgebroid structure over 
A= 4>~ l {l} = R x TM 2 x R 3 . Note that the Lie affgebroid structure oni = Mx TK 2 x R 3 is a 
special type of Lie affgebroid structure called Atiyah affgebroid structure (see Section 9.3.1 in 
for a general construction). Moreover, the affine bundle : Ix TM 2 xl 3 — > R 3 is modelled 
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on the vector bundle t v : V = _1 {O} = K x TR 2 x R 3 -> R 3 . Thus, (t, a;, y; i, y, w x , w v , w a ) 
may be considered as local coordinates on .4 and V. 
In this case, from ([6]), we have that a basis of sections of T A A is 
r - <9~ d ~ d 

| T = (e , — ), Ti = (ei, — ), T 2 = (e 2 , ^), T 3 = (e 3 ,0), T 4 = (e 4 ,0), T 5 = (e 5 ,0), 

Consider now the following mechanical problem. A (homogeneous) sphere of radius r > 0, mass 
m and inertia about any axis k 2 , rolls without sliding on a horizontal table which rotates with 
time-dependent angular velocity £l(t) about a vertical axis through one of its points. Apart 
from the constant gravitational force, no other external forces are assumed to act on the sphere. 
Therefore, the Lagrangian of the system corresponds to the kinetic energy. Moreover, observe 
that the kinetic energy may be expressed as a Lagrangian L : A — > M: 

L(t,x,y;x,y,uj x ,ujy,uj z ) = ^(mx 2 + my 2 + k 2 (u 2 x + u; 2 +w 2 )), 

where (lo x , u> y , ui z ) are the components of the angular velocity of the sphere (see 0], for more 
details). 

After some straightforward calculations using (|10p , we deduce that the Poincare-Cartan sections 
associated with L are given by: 

Q L = -Lfa + mxf 1 +myf 2 + k 2 uj x f 3 + k 2 uj y f 4 + k 2 uj z f 5 , 

Q L = <f AA L A <£o + mf 1 A V 1 + mf 2 A V 2 + k 2 f 3 A V 3 + k 2 f 4 A V 4 + k 2 f b A V 5 
+k 2 u x f 5 AT 4 + k 2 uj y f 3 AT 5 + k 2 uj z f 4 A f 3 . 

Since the ball is rolling without sliding on a rotating table then the system is subjected to the 
affine constraints: 

'J 1 = fi(i)y + x- rw y , 
* 2 = -Q(t)x + y + ruj x , 
which define an affine subbundle B of A. Then, we have that 

d T A A^,i = ci'(t)y^ +n(t)f 2 + V 1 ~rV 4 , 

d T A A^2 = _ fi '( f ) a; o _ 1 + V 2 + r y3 

Thus, the Lagrange-d'Alembert equations (see ([55]) ) for the system determined by (L,B) are: 



•mx — 


-Ai, 


my = 


-A 2 , 


< k 2 LU x — 


-rA 2 


k 2 Uly = 


rAi, 


k 2 u z = 


o, 



and the constraints = and \& 2 = 0. 

In this basis of sections of vector subbundle F C T a A\b — > B is given by 

r 1 ~ r ~ 1- r ~ 1 

{ Zi = 7l + y^y 4 , ^2 = V 2 - j=V a \ 

v m k z Tfi k ) 
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and then, the matrix C is 



The projector P over T A B is 



P = Id + 



mk 2 



( _l_r!. 

m k 2 




;Z 1 ® d^^ 1 + 




mk 



k 2 + mr 

Since the solution of the unconstrained dynamics is 

R L =f + xf X + yf 2 + LU x f 3 + LUyfi + LJ z f 5 , 

then the solution of the nonholonomic problem will be 

Rnh = P(Rl) =T + xTi + yT 2 + uj x T 3 + io y T± + lu z T 5 
+ 



k 2 + 



TY1 TD A 1 ^ 

' ' (fi'(i)j/ + Cl(t)y)Zi - —t- -(n'(t)x + n(t)x)Z 2 . 



mr* 



k 2 + 



mr* 



Thus, 



d 



a 



d 



P T f{Rnh) = ^+x— + y— 



dt dx 
mr 

+ 



k 2 + mr 2 



dy k 2 + mr 

(n'(t)x + n(t)x) 



d 



+ 



duj x k 2 + mr 2 



dx k 2 + mr 

TT) v n 

' (fl'(t)y + n(t)y) 



''" (fi'(i)z + fi(t)±)^ 



duj„. 



and, therefore, the equations of motion of the nonholonomic problem are 

k 2 



k 2 



-(n'(t)y + n(t)y), 

{n'(t)x + Q(t)±), 

(n'(t)y + Cl(t)y), 



Now, we take the coordinates (t,x,y;pt,p x ,p y ,u x ,u y ,u z ) on A + and the corresponding coor- 
dinates (t,x,y;p x ,p y ,u x ,Uy,u z ) on V* . Then, we obtain that the extended Legendre transfor- 
mation LegL ■ A — > A + and the Legendre transformation legL '■ A — > V* associated with L are 
given by 

Leg L (t,x,y;x,y,u> x ,u> y ,uj z ) = (t, x, y; —L, mx, my, k 2 uj Xl k 2 u> yi k 2 u z ), 
leg L (t,x,y;x,y,uj x ,uj y ,uj z ) = (t,x,y;mx,my,k 2 uj x ,k 2 uj y ,k 2 uj z ). 





k 2 


y 


k 2 + mr 2 


< . 


mr 


u x 


k 2 + mr 2 




mr 


LOy 


k 2 + mr 2 


k w 2 


= 0. 
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Thus, it is easy to see that the image of B under the Legendre transformation, which will be 
denoted by B, is defined by the vanishing of the functions 

1 r 
ip 1 = n(t)y + —p x - —u y , 



m 



ijj 2 = -Q.{t)x H p y + -r^u x . 

m k z 

Moreover, the induced Hamiltonian section h — Leg^ ° leg^ 1 : V* — ► A + is given by 

h(t,x,y;p x ,p y ,u x ,u y ,u z ) = (t,x,y; -H,p x ,p y ,u x ,u y ,u z ), 
where the function H is 

H(t,x,y;p x ,p y ,u x ,u y ,u z ) = \ {^(pl+Py) + ^( u l + u l + u l)) ■ 

A curve 7 : I — > V* , t (i, x(i), y(t);p x (t),p y (t), u x (t), u y (t), u z (t)), is a solution of the 
equations of motion for the nonholonomic Hamiltonian system (h, B) if it satisfies the equations 



Px 


= -Ai, 


Py 


= -A2, 


u x 


= -r\ 2 


iiy 


= r-Ai, 




= 0, 



and the constraints ip 1 = and t/j 2 = 0. 

In this particular example, the expression of the nonholonomic bracket is given by 



{h',h"} nh 



d(H'-H") I QH' dH" _ dH' dH" \ , / dH' dH" _ dH' dH" 



dt 



\ dx 



dx 



\ dy dp y 



, f dH'dH" dH' dH" \ ( dH' dH" dH' dH" \ f dH' dH" dH'dH" \ 

' ' 1 : ' '- ) vydux du z du z du x ) z \du y du x du x du y ) 



'' y du z du y du y du z 



k m 



k 2 +r 2 m 



k m 



xdH" 



dH" 



dH"' 



dp 



d(H'-H) _ r d(H'-H) 
dp x du y 



k 2 -\-r 2 m 



Q'(t)x fi(i)^ m d Q y k 2 ^ u v 9 dZ z u z d d H u v > 



d(H'-H) | r d(H'-H) 



dp y 



du x 



+ 



k m 
k 2 +r 2 m 



W)v - k% + + - u x mL) 



d(H"-H) d(H"-H) 
dp x du v 



k m 
k 2 -\-r 2 m 



d(H"-H) d(H"-H) 
dp y du x 



m dx ~ r "Wa p „ ^W\ az du x du 

v( t )x-m^-^-M^-^) 

for sections h! and h" with associated functions H' and H", respectively. 

VII. Conclusions and Future work 

We have developed a general geometrical setting for nonholonomic mechanical systems in the 
context of Lie affgebroids. We list the main results obtained in this paper: 

• The notion of regularity of a nonholonomic mechanical system with affine constraints 
on a Lie affgebroid was elucidated and characterized in geometrical terms. 
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• In the regular case, the solution of the nonholonomic problem was obtained by pro- 
jecting the unconstrained one using several decompositions of the prolongation of the 
Lie affgebroid along the afhne constraint subbundlc. 

• The hamiltonian formalism for nonholonomic systems is completely analyzed and a 
nonholonomic bracket was defined. 

• Several examples were discussed showing the versatility of this geometric framework. 

In a forthcoming paper we will study the reduction of the Lie affgebroid nonholonomic dyna- 
mics under symmetry and we will obtain a Lie affgebroid version of the momentum equation 
introduced in [7] for Lie algebroids. 

Other goal we have proposed is to develop a geometric formalism for vakonomic Mechanics and 
optimal control theory on Lie affgebroids. 

In addition, we will explore the construction of geometric integrators for mechanical systems 
on Lie affgebroids and, in particular, for nonholonomic systems. 
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